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The  main  objective  of  this  study  was  to  explore  a  new 
technique  to  improve  the  dynamic  response  and  vibration 
characteristics  of  composite  plates  or  beams  by  using  shape 
memory  alloy   (SMA) .  Detailed  descriptions  of  the  making, 
testing,   and  modeling  of  SMA  hybrid  structures  are  reported 
in  this  paper. 

The  unique  properties  of  SMA  and  its  applications  are 
described  at  the  beginning.  Analyses  were  made  to  determine 
the  transient  thermal  response  of  SMA  laminates  activated  by 
electric  current.  A  finite  element  analysis  to  verify  the 
activation  of  beam  specimens  is  presented  in  the  paper  and 
shows  good  agreement  with  the  test.  A  complete  tooling, 
fabricating,  and  wire-straining  device  of  the  specimens  is 
also  developed. 
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An  analytical  tool  was  developed  to  calculate  such 
properties  of  SMA  hybrid  laminates  as  the  A,  B,  D  matrices 
of  in-plane,   coupling,  and  bending  stiffness. 

A  theoretical  method  of  predicting  the  internal  stress 
distribution  of  SMA  hybrid  laminates  was  also  formulated 
under  a  fully  activated  condition  of  SMA  wires.  The 
classical  lamination  theory  (CLT) ,   combined  with  this 
method,  gives  a  feasible  approach  to  estimate  stress-strain 
relations  in  the  SMA  hybrid  laminates. 

The  dynamic  responses  of  SMA  hybrid  laminated  beams 
(thin  laminates)  were  tested  experimentally.  The  results 
showed  evidence  of  natural  frequency  shifting  for  beams  with 
clamped- clamped  boundary  condition.  The  failure  rate  of  the 
experiment  was  high  due  to  the  particular  boundary 
condition,  and  other  boundary  conditions  were  not 
investigated. 

Finally,  a  dynamic  model  based  on  large  deformation 
theory  was  derived  for  the  SMA  hybrid  composite  sandwich 
beam.  This  model  can  predict  the  natural  frequencies  at 
first  order  nonlinearity  and  the  dynamic  response  under 
harmonic  excitations.  Analytical  results  showed  good 
improvements  in  the  beam  vibration  suppression  under  hinged- 
simply  supported  boundary  condition.   Parametric  studies 
revealed  that  the  in-plane  forces  induced  by  embedded  SMA 
wires  in  composites  were  the  major  contributor  in  beam 
vibration  control  under  the  said  boundary  condition. 
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CHAPTER  1 
INTRODUCTION 


1.1    Backgrounds  of  The  Project 

There  have  been  many  problems  in  dealing  with  vibration 
control  in  every  engineering  field.  The  relevant  techniques 
in  solving  vibration  problems  fall  into  two  major  groups; 
passive  vibration  control  and  active  vibration  control. 
Passive  vibration  control  has  been  widely  used  in  structure 
components  of  spacecraft,   aircraft,   ships,  ground  vehicles, 
rotating  machines,  buildings,  bridges,   and  so  on.  Simply, 
passive  vibration  control  uses  tuned  masses,  viscous 
dampers,   dashpots,  or  visco-elastic  tapes  fixed  on  the 
structure's  components  to  alter  the  damping  level  and  the 
system  response.  Active  vibration  control,  on  the  other 
hand,  uses  sensors  to  measure  the  structure's  response  and 
generate  signals  to  the  feed-back  control  block  to  activate 
actuators   [1] .  Materials  having  both  sensors  and  actuators 
in  closed  loop  control  are  called  adaptive  materials.  The 
active  vibration  control  systems  possess  both  sensors  and 
actuators  which,   fastened  on  structure  components, 
constitute  what  is  called  an  adaptive  structure. 

New  kinds  of  materials,   called  smart  materials,  have 
come  to  play  a  more  important  role  in  the  last  decade.  Smart 
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materials   [1]   have  unique  abilities  to  change  their  physical 
properties  and  geometries   (i.e.,   stiffness,  damping, 
viscosity,   shape,  etc.)   in  response  to  certain  stimuli.  At 
this  time,  optical  fibers,  piezoelectric  polymers, 
piezoceramics,  electro-strictors,  electro-rheological 
fluids,   shape  memory  alloy  (SMA) ,  and  magnetostrictors  are 
examples  of  some  smart  materials  most  often  studied.   It  is 
said  that  adaptive  materials  and  structures  are  a  subset  of 
smart  materials  and  structures   [1] . 

Smart  structures   [1]   have  embedded  or  bonded  intrinsic 
sensors  which  recognize  and  measure  the  intensity  of  the 
stimuli.  These  stimuli  may  be  stresses,  strains, 
temperatures,  electrical  voltages,  magnetic  fields,  chemical 
compounds,  or  radiation.  Smart  structures  have  embedded  or 
bonded  intrinsic  actuators  that  can  respond  to  stimuli. 
Smart  structures  have  certain  mechanisms  for  controlling  the 
responses  to  the  stimuli  according  to  a  predetermined 
relationship  and  are  capable  of  selecting  the  best  response 
if    more  than  one  option  is  available. 

In  the  last  decade,  adaptive  structures  have  shown 
great  potential  for  modal  control  of  large  space  struc-  - 
tures  or  for  use  in  certain  mechanical  devices.  Smart 
materials  are  employed  to  reach  full  adaptive  control  of  the 
dynamic  responses  of  structures,   such  as  space  cranes,  robot 
arms,   and  similar  mechanisms.  Smart  materials  have  not  only 
reduced  the  weight  of  traditional  adaptive  actuators  but 
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have  stimulated  the  evolution  of  a  new  concept  of  dynamic 
property  tuning  in  smart  structures . 

The  shape  memory  alloy,  used  in  the  past  decade  on 
robot  arms   [2] ,  has  proven  to  be  an  effective  actuator 
because  of  its  simple  construction  and  smooth  mobility.  The 
utilization  of  SMA  enables  the  designer  to  reduce  the 
complexity  of  mechanical  elements.   In  addition,   the  facts 
that  it  has  long  fatigue  life  and  corrosion-resistant 
capability  have  made  the  material  more  useful  to  engineering 
applications.  To  activate  the  SMA  requires  heat,   a  common 
form  of  energy.  Heat  can  be  generated  by  electric  current 
through  SMA  wires  in  the  same  way  as  electric  current 
through  electric  resistors  generates  heat.  Through  thermal 
activation,   the  SMA  hybrid  structures  can  change  its 
stiffness  and  material  strength  as  well  as  the  structure's 
natural  frequencies. 

A  large  structure  has  many  spars,   rods,   and  beams.  For 
instance,   the  structure  may  be  a  truss  or  a  drawbridge.  The 
modal  response  of  a  drawbridge  would  change  if  only  one 
property  (e.g.,   stiffness,   strength)  of  one  element  of  the 
structure  changed.  A  light  structure  possessing  many  members 
may  not  be  strong  or  stiff  enough  to  withstand  changing 
external  loads,  or  a  stiff  structure  is  strong  but  too 
heavy.  By  replacing  a  few  of  the  components  with  SMA  hybrid 
beams,  one  could  achieve  the  purpose  of  improving  the 
dynamic  responses  of  the  structure.  Therefore,  because  of 
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weight  saving,   it  could  be  more  efficient  to  use  the  SMA 
material  for  dynamic  control . 

1.2     Special  Features  of  SMA 

There  are  two  types  of  solid  state  transformation  of 
materials.  One  is  diffusional  transformation,  where 
migrating  atoms  randomly  form  a  new  phase  with  a  different 
chemical  composition  from  the  original  matrix.  This 
progression  of  transformation  is  dependent  upon  both  time 
and  temperature.  The  other  transformation  is  displacive 
transformation .  which  does  not  require  either  atomic 
migration  or  change  in  the  chemical  nature  of  the  matrix  to 
alter  its  solid  state.   Instead,   its  crystal  structure  is 
merely  rearranged  into  a  more  stable  one.  This  displacive 
progression  of  transformation  is  only  dependent  upon 
temperature    [2] . 

An  example  of  displacive  transformation  is  the 
austenite  ^  martensite  phase  transformation  of  SMA.   The  2-D 
diagram.   Figure  1.1  schematically  shows  the  progression  of 
transformation  from  one  phase  to  the  other.  The  volume 
fraction  of  martensite  increases  as  temperature  decreases  at 
the  parent,   or  austenite  phase.   In  the  austenite  phase,  each 
atomic  layer  is  packed  together  in  a  rectangular  shape,  and 
each  fine  layer  tilts  to  achieve  the  twinned  martensite 
phase.  This  phenomenon  is  called  Bain  strain   [2] . 
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Figure  1.1.        Displacive  transformation  of  austenite  to 
martensite  [2] 

Nickel -titanium  (Ni-Ti)   alloy  (Nitinol)   is  a  kind  of 
polycrystal  SMA.  The  metal  crystal  structure  has  grain 
boundaries  in  both  phases   [3] .   It  is  stiff er  than  most  of 
the  SMA  products  currently  available.  This  particular 
material  is  applicable  in  the  fields  of  actuators, 
electrical  connectors,  and  fasteners.   It  also  has  superior 
characteristics  in  corrosion-resistance. 

In  general,   this  Nitinol -type  of  shape  memory  alloy 
will  undergo  an  accommodation  step.   The  change  in  lattice 
shape  is  not  merely  caused  by  Bain  strain  because  the  other 


layers  of  austenite  may  not  always  change  to  the  slant 
typical  of  martensite,   shown  in  Figure  1.1,  but  may 
accommodate  themselves  to  a  more  relaxed  shape,  dubbed 
twinning .  Figure  1.2  shows  the  general  shape  of  martesite 
after  cooling,  and  the  shape -changing  process  is  reversible 
This  phenomenon  is  called  lattice  invariant  shear   [2] . 


Figure  1.2.        Twinning  of  shape  memory  alloy 


The  twinning  process  plays  a  key  role  in  the  shape 
memory  effect.  We  can  see  from  Figure  1.2  that  the  boundary 
at  each  fine  plate   (four  thin  plates  above  and  four  thin 
plates  below  the  midplane   [2] )   slants  at  a  different  angle, 
but  still  has  all  the  atomic  bonds  intact  without  slippage. 
The  energy  state  of  this  condition  is  very  low,   and  it 
becomes  very  mobile. 
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An  optical  micrograph  by  Duerig  et  al .    [2]  of  surface 
relief  due  to  the  formation  of  martensite  in  a  Cu-Zn  shape 
memory  alloy  shows  the  twinning  stage    [2]   and  better 
depicts  the  true  twinning  shape  in  contrast  to  the  2-D 
lattice  diagram  shown  previously  in  Figure  1.2.  This 
phenomenon,  true  in  fine  plates  of  molecular  structures, 
is  also  true  in  thicker  plate  sizes  up  to  the  entire  crystal 
structure . 

Twinning  in  ordered  alloys  is  not  always  impossible,  as 
Cu-Al-Ni  martensite  shows  both  {l2l}^,  and  {lOl},,  twinning 
[3] .   Cu-14Al-3Ni  is  a  type  of  single  crystal  SMA,   where  the 
stress-assisted  transformation   (stress  induced  martensitic 
transformation)  produces  a  typical  example  of  memory  effect. 
At  T<Mf     (where  T  is  the  temperature  at  which  stress  is 
applied,        is  martensite  final  temperature) ,   the  true 
activation  temperature  of  130°C  is  higher  than  that  of  non- 
stress  induced  martensitic  transformation.  Yet,   at  M3<T<Af 
(where       is  the  martensite  starting  temperature  and  A^  is 
the  austenite  final  temperature) ,   the  true  activation 
temperature  becomes  lower   (approximately  100°C)  .  Therefore, 
the  twinning  can  be  easily  changed  into  deformed  martensite 
alloy  by  applying  stresses  or  by  heating  up  to  austenite 
alloy  as  shown  in  Figure  1.3.   For  example,   a  twinned 
martensite  SMA  wire,   like  Cu-Zn  single  crystal  material, 
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Figure  1.3.     Shape  memory  process  shown  microscopically 

deforms  plastically  (pseudoelastic  deformation)   under  a 
tension  load  above        (martensite-starting  temperature) .  The 
deformed  crystal  length  tends  to  shrink  when  heat  is 
applied,  as  if  the  heat  energy  is  absorbed  by  the  crystal  to 
pack  itself  more  tightly  to  achieve  a  higher  energy  state 
(austenite  state) .  This  phenomenon  is  the  memory  effect 
found  in  most  of  the  SMA  materials.   In  addition,   the  single 
crystal  SMA  usually  has  a  better  memory  effect  than  that  of 
a  polycrystal  SMA   [2]  , 

1 

1 

This  phenomenon  clearly  explains  why  Nitinol  wire, 
which  deforms  plastically  above  M3   (approximately  70 °F)  at 
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a  constant  load  "P, "  must  shrink  back  to  austenite  phase 
when  heat  is  applied.   It  also  indicates  that  properties  of 
the  wire  will  have  changed  at  the  activation  temperature. 

1.3    Outline  of  This  Study 

The  research  [1]   reports  that  SMA  wires  can  be  embedded 
in  various  kinds  of  composite  materials,   resulting  in 
vibration  suppression,   to  change  structure  stiffness  to 
avoid  resonance   [1,4],   for  acoustic  vibration  control  [5], 
and  for  working  as  a  geometry  variation  device. 

The  objective  of  the  study  is  to  seek  an  effective  way 
to  suppress  vibration  by  using  SMA  as  an  actuator  in  a 
composite  structure.  Five  major  areas  of  studying  SMA  hybrid 
composites  are  described  in  this  paper.  These  are:  [Thermal 
response  of  SMA  hybrid  composites;  The  fabrication  of  SMA 
hybrid  specimens;  Stress -strain  in  SMA  hybrid  laminate; 
Dynamic  test  of  SMA  hybrid  laminates;  The  dynamic  model  of 
SMA  sandwich  structure] . 

^•3-1 — Thermal  Response  of  SMA  Hybrid  Composites 

Analytical  means  were  introduced  [6]   by  employing  the 
Finite  Difference  Method  to  analyze  the  transient  thermal 
response  of  a  SMA-graphite/epoxy  beam.  The  beam  temperature 
is  the  only  key  related  to  the  activation/deactivation  of 
the  hybrid  structure.  The  author  uses  the  thermal  analysis 
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module  in  ANSYSS  based  on  the  modeling  method  described  by 
Hawkins  et  al .    [6],   to  investigate  the  temperature  in  the 
beam  cross  section.  The  specimens,  made  of  SMA-fiber 
glass/epoxy,  were  also  tested  to  compare  with  the  analytical 
thermal  response.  Results  are  reported  in  Chapters  2  and  4. 

1 • 3 ■  2  The  Fabrication  of  SMA  Hybrid  Specimens 

The  SMA  wire  embedding  concept  is  dscribed  by  Rogers 
and  Barker   [4] .  Based  on  this  concept,   the  author  developed 
a  different  wire-straining  and  fabricating  method  for  the 
SMA  hybrid  laminate.  Details  is  presented  in  Chapter  3. 

1-3.3  Stress-Strain  in  SMA  Hybrid  Laminatp 

Because  the  hybrid  laminate  performs  differently  from 
that  of  the  general  composite  laminate,   the  hybrid  lamina 
properties  must  be  calculated  in  a  different  fashion   [7] . 
The  compound  effect  of  the  composite  thermal  stress,  the 
boundary  stress,   and  the  special  recovery  stress    [2]   of  SMA 
wires  produce  the  laminate  in-plane  force   [7]  when  the 
activation  occurrs.  As  a  result,   the  combined  effects  of 
both  in-plane  force  and  internal  stress  at  the  SMA/composite 
interface  present  a  complicated  stress-strain  problem.  The 
author  made  a  first  attempt  to  estimate  the  internal  stress 
distribution  by  using  the  principles  of  solid  mechanics  as 
Detailed  in  Chapter  5. 


ANSYS  is  a  finite  element  package  for  structural,  and 
thermal  analysis 
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1-3 -4  Dynamic  Test  of  SMA  Hybrid  Laminates 

The  test  work  by  Rogers  and  Barker  [4]  and  Saunders  et 
al.    [5]   showed  the  phenomenon  of  natural  frequency  shifting 
at  the  activation  temperature.  The  author  uses  the  same 
concept,  with  different  SMA-composite  laminate,   to  verify 
the  natural  frequency  change.  The  modal  test  results  are 
listed  in  Chapter  6. 

1-3 -5 — The  Dynamic  Model   of  SMA  Sandwich  Structure 

There  are  several  studies  of  SMA  hybrid  laminated  beams 
and  plates,  but  none  of  them  has  involved  the  employment  of 
SMA  in  a  sandwich  structure.  The  author  used  the  nonlinear 
approach,  based  on  the  proven  theory,   to  find  the  governing 
equations  of  the  SMA  composite  sandwich  panel.  The  complete 
derivations  and  the  analytical  results  are  discussed  in 
Chapter  7. 


CHAPTER  2 
SMA  HYBRID  COMPOSITE  PLATE 

2.1    Application  of  SMA  Hybrid  Laminates 

Vibration  suppression  can  be  achieved  by  using  the  SMA 
active  strain  energy  tuning   (ASET)   method   [1,4],   where  the 
SMA  wires  are  plastically  elongated   [2]   prior  to  being 
embedded  in  a  composite  structure.  This  gives  the  structure 
the  ability  to  create  large  distributed  loads  in  the 
composite  material  upon  heating.  The  wires  are  generally 
placed  in  the  structure  in  such  a  way  that  no  bending 
deflection  will  occur  during  the  activation  (heating) 
process.   Instead,  the  structure  is  placed  in  a  residual 
state  of  strain.  The  resulting  stored  energy  changes  the 
modal  response  of  a  structure  in  a  way  similar  to  that  of  a 
guitar  string   [4] .  The  result  is,   the  natural  frequencies  of 
the  hybrid  laminate  increase  at  the  activation  temperature. 
This  is  a  direct  use  of  shape  memory  effect  as  described  in 
Figure  1.3  of  Chapter  1 . 

2.2  ASET  Set -Up  Concept 

The  author  used  a  specimen   [0/90/SMA/0/90] 3  SMA-GF/EP 
hybrid  composite  laminate  beam   (GF/EP:  glass  fiber 
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unidirectional  prepreg) .  The  embedded  SMA  was  heated  up  to 
accomplish  the  transformation  from  the  martensitic  phase  to 
the  austenite  phase.  During  this  procedure,   the  Young's 
modulus  changes  by  a  factor  of  three.  The  heating  process  is 
accomplished  by  running  an  electric  current  though  the  wire. 
This  procedure  is  shown  schematically  in  Figure  2.1. 

The  consolidation  process  of  composite  laminate  in 
fabrication   [4]   can  affect  the  mechanical  behavior  of  the 
SMA  actuators,  because  the  wires  are  constrained  from  moving 
during  the  cure  cycle  by  clamping  each  wire  to  the  curing 
plate.  A  clamped- clamped  hybrid  composite  beam   (Figure  2.2), 
has  a  constant  current  power  supply  controlled  manually  to 
bring  the  temperature . of  the  beam  (activation)   to  the 
desired  austenite  finish  temperature. 

An  electromagnetic  displacement  transducer  is  used  to 
sense  the  dynamic  response  of  a  composite  beam.   It  can 
convert  analog  data  to  digital  data  and  transmit  it  to  a 
data  acquisition  system  for  modal  testing   [4,5]  . 

2.3     Thermal  Response  of  SMA  Hybrid  Beam 

The  author  chose  Nitinol  wire  as  the  SMA  wire  for  use 
in  the  hybrid  composites.   The  SMA  metal  content  is  55% 
nickel  and  45%  titanium.   The  SMA  density  is  6.45  g/cm\  The 
Melting  point  of  SMA  is  1240°C-1310°C.  The  SMA  Young's 
modulus  Eg  for  the  austenite   (parent  phase)    is  11.0E06  psi. 
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The  SMA  Young's  modulus  E„  for  martensite   (twined  material 
phase)   is  3.5E06  psi .   The  SMA  Wire  diameter  is  0.01".  Data 
source  is  from  the  manufacturer,   Precision  Wire  Corporation 
[8]  . 

The  SMA  wires  are  heated  to  activate  and  cooled  to 
deactivate.  The  speed  of  adaptability  is  governed  solely,  by 
the  thermal  activation  and  deactivation  times  for  the  SMA 
hybrid  composite  beam.  Hawkins  et  al .    [6]   reported  that  an 
impulse  electric  power  supply  can  effectively  activate  a 
hybrid  beam  in  a  short  time.  The  energy  supply  can  then  be 
cut  back  to  sustain  the  activation  temperature. 
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Because  of  the  strong  relations  between  the  activation 
time  and  the  dynamic  response,  a  thermal  analysis  of  the 
hybrid  composite  is  necessary  [6] .  Furthermore,   a  series  of 
tests  on  real  specimens  is  imperative  to  verify  the 
transient  thermal  response  inside  the  hybrid  composite  beam. 

The  theoretical  approach  to  transient  thermal  response 
[6]   is  shown  in  Figure  2.2.  Equation  (2.1)  describes  heat 
conduction  in  the  hybrid  laminate,  with  heat  source  q 
(generated  by  electric  resistance  heating)   of  a  SMA  wire. 
The  g  (volumetric  heat)   is  the  energy  that  functions  as  an 
active  force  for  the  SMA  actuator. 

^^"-^  ^  -5^  ^^^-3y '     -^'^-Tt  (2.1) 

Equation   (2.2),  where 
I ( t )   =    current   ( amp ) ; 

Re  =  electric  resistance   (ohm)  ,  p^*h  /  1^ ; 

Pe  =  electric  resistivity; 

V„  =  volume  of      wire   (mM   =  l^  L; 

1^  =  H  SMA  cross  sectional  area; 

L     =  SMA  unit  length   (set  to  1"); 

y  ~   t7         -   '       (the  actuator)  (2.2) 

=  0,  (the  composite) . 
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Figure  2.2.  Cross  section  of  a  hybrid  laminate  and  a  thermal  model 
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Choosing  the  following. non-dimensional  terms,  we  obtain 


range  from  0  to  1; 
range  from  0  to  a; 
range  from  0  to  b; 

real  time   (hr.  or  sec . ) ; 

non-dimensional  time; 
Martensite  final  temperature; 
Austenite  final  temperature; 
ambient  temperature. 


(2.3) 


e 


(2.4) 


'  =    "       V  -  y 


(2.5) 
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K=  K>   thermal  conductivity  in  actuator   (SMA)   at  0°, 

=  ^cxi   thermal  conductivity  in  composite  at  0°, 
ky  =  k„,   thermal  conductivity  in  actuator  (SMA)   at  90°, 
ky  =  key,   thermal  conductivity  in  composite  at  90°, 
Cp  =  Cp„,  heat  capacity  in  actuator   (SMA)  , 
=  Cpc,  heat  capacity  in  composite. 


By  substituting  the  previous  terms  in  equation  (2.1),  the 
equation  is  converted  to  the  non-dimensional  equation 


the  boundary  conditions  and  initial  conditions,  according  to 
Figure  2.2,  are 


=  0  ,  at  X=0,  A,  where  A=^. 


and 


(2.9) 


=  0  '  at  y=0  (2.10) 


^  =  B^e  ,  at  y=B,  where  B=^  (2.11) 
6=0,  at  T  =  0  (2.12) 


-  =  Biot  number.  (2  13) 
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The  Biot  number  represents  the  heat  convection  from  the 
beam  surface  to  the  ambient  air  {T„  =  70 °F)  .  The  non- 
dimensional  equation  can  be  solved  by  the  Finite  Difference 
Method  or  by  a  commercially  available  thermal  analysis 
package.  The  author  used  ANSYS  to  verify  the  transient 
thermal  response  for  the  experimental  specimen. 

For  comparison  purposes  in  thermal  analysis,  a 
verification  experiment  measuring  the  transient  thermal 
response  of  using  an  experimental  specimen  is  very 
important.  Analysis  of  the  said  heat  conduction  equation 
(2.8)   can  explicitly  show  us  the  nondimensional  temperature 
rise  and  fall  during  SMA  wire  activation-deactivation 
process  Hawkins  et  al . [6] .  The  solution  can  also  show  that 
the  required  electric  power  can  rapidly  activate  the  beam, 
and  the  vibration  control  may  be  achieved  in  a  proper  time 
span.  Theoretically,   an  electric  impulse  can  quickly 
activate  the  hybrid  beam  ( [6]   at  various  Biot  numbers)  and 
at  different  beam  constructions.  However,   the  technique  for 
deactivating  a  hybrid  structure  may  pose  a  problem,  because 
the  ambient  environment  is  potentially  changeable  that  a 
cost-effective  heat  dissipation  device  can  be  difficult  to 
design. 


CHAPTER  3 

THE  SPECIMEN  FABRICATION  AND  TESTING  PLAN 
3.1     Specimen  Fabrication  of  SMA  Hybrid  Beam 

The  reason  behind  the  use  of  SMA  wires  in  composite 
laminates  is  its  ability  to  increase  stiffness  in  composite 
beams.  The  author  used  composite  laminates  because  the  SMA 
hybrid  aluminum  plate  or  hybrid  steel  plate,  made  with 
aluminum  or  steel  mixed  with  the  SMA  wire,   has  too  high  a 
melting  point  and  the  properties  of  the  wire  are  completely 
destroyed,  because  SMA  no  longer  preserves  its  quality  at 
the  martensite  phase  or  remembers  its  features  at  the 
austenite  phase.   The  curing  temperature  for  the  graphite/ 
epoxy  composite  is  only  350°F   (177°C)  ,   so  the  difference 
between  room  temperature  and  curing  temperature  (autoclave) 
is  only  about  250°F  to  300°F   {120°C  to  150°C)    [3]  .  Therefore, 
they  are  still  on  the  safe  side  of  SMA  operating  range, 
which  is  around  170°C,   and  the  properties  will  scarcely  be 
affected  in  the  fabrication  process. 

In  the  experiments,   we  found  some  uncertainties  in  the 
properties  of  the  hybrid  beam  before  and  after  fabrication. 
For  the  first  batch  of  specimens,   the  wire  was  fixed  on  the 
tool  plate   [6]   and  not  strained  to  5-6%  as  recommended  by 
Rogers  et  al .    [7] .  Trial  and  error  taught  us  not  to 
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overstrain  the  wires,   instead  to  permit  them  to  resume  the 
austenite  phase  first  and  cool  before  embedding  in  prepregs. 
The  second  batch  of  specimens  were  made  by  a  modified  wire- 
straining  tool  and  fabrication  tool,  which  can  strain  the 
wires  to  10%,  and  this  was  found  to  be  satisfactory. 

3.1.1     Selection  of  Fabrication  Tool 

The  tool  materials  and  how  they  are  set  up  are  very 
important  factors  affecting  the  produced  specimen.  First, 
the  author  listed  some  important  physical  properties  and 
mechanical  properties  of  the  Nitinol  SMA.  The  following  data 
are  provided  by  the  Precision  Wire  Company   [8] : 


3.1.1.1     Properties  of  Nitinol 
Physical  properties 

Density  

Electrical  resistivity  

Coefficient  of   

Thermal  expansion   

Thermal  Conductivity  

Mechanical  properties 

Ultimate  Tensile  Strength.. 
Yield  Strength   


6.45  g/cm^    0.234  Ib/in^ 
76  microhm-cm  (Martensite) 
82  microhm-cm  (Austenite) 
6.6xlO-V°C  (Martensite) 
11.0xlO-V°C  (Austenite) 
0 . 17  watt/cm- °C     (K^  ) 

125  ksi 

80  ksi  (Austenite) 
10  ksi  (Martensite) 
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Elongation   60  percent  (Martensite) 

20  percent  (Austenite) 

Young's  Modulus   3.5x10*  psi  (Martensite) 

11.0x10*  psi  (Austenite) 

Hardness   42  HRA 

Poisson's  Ratio   0.33 

Other  factors 

Stock  size  used  (in  shop)..  0.01"  diameter  wire 

(micrometer  measured) 

Cost/per  foot   approximately  $1.31  (1992) 

Martensite  final   (  Mf  )  . . . .  approximately  20°C 

Austenite  start     (  A^  )  . . . .   38°C  (100°F) 

Austenite  final     (  A^  )  .  .  .  .   130°F  to  290°F 
Thermal  Expansion  Coefficient  of  GF/Ep 

0(1= .  045xlO'V°C 

Qf2=20 .  2xlO"V°C 

The  hybrid  laminate  specimen   [0/90/SMA/0/90] 3,   shown  in 
Figure  3.1,   is  2.87cm  x  2ncm  beam  (first  batch)   for  the 
transient  thermal  test  and  the  dynamic  test. 

The  base  tool  plate  and  the  Nitinol  wire-straining 
frame   (Figure  3.2)   should  be  steel,  because  steel  has 
approximately  the  same  thermal  expansion  coefficient  as  that 
of  the  austenite  SMA.   Therefore,   steel  alloy  with  thermal 
expansion  coefficient  of  llxlO-V°C  and  aluminum  alloy  with 
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thermal  expansion  coefficient  of  22xlO"V°C)  were  used  as  the 
tooling  materials. 


cul    at  1    I  Ik'  I  ot  • 


Can    be    made    4    speci  men  r" 

A-  A 


Figure  3.1.   Specimen  fabrication  diagram   (type  1  - 
2,3;   type  2  -  no .   1 ) 
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3.1.2     Tool  Design  and  Fabrication  Procedure 
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The  author  used  a  1.25  X  1.25  X  O.l"  aluminum  extruded 
L-shape  beam,  cutting  teeth  on  it  with  a  rotary  hacksaw.  The 
teeth  are  used  for  aligning  SMA  wires  preparatory  to  burying 
the  wires  in  the  composite  beam.  Grinding  and  refining  the 
teeth  were  done  with  a  DREMEL  miniature  tool    (large  and 
small  collet  with  various  kinds  of  engraving  bits  or  sanding 
bits  to  trim  the  teeth  on  the  L-shape  beam  in  Figure  3.2) . 

Two  aluminum  dams  were  cut  to  shape,   and  each  was 
placed  at  the  corner  of  the  L-shape  beam,   forming  a 
rectangular  frame,  as  seen  in  Figure  3.2.  Then,  we  aligned 
the  four  corners  to  drill  through  the  L-shape  beam  and  dams. 
Tapped  the  dams  with  internal  threads  for  screw  locking  to 
strain  the  SMA  wires  before  fabrication.  Secondary  base 
plates  were  cut  and  bonded  to  the  right  hand  side  of  dams  to 
fill  the  gap  between  the  dams  and  the  base  plate. 

The  secondary  base  plate  was  trimmed  to  shape  and  put 
on  the  base  plate  to  complete  the  basic  frame  assembly. 
Figures  3.2  and  3.3  show  the  perspective  view  and  cross - 
sectional  view.  Each  end  of  the  dam  was  chamfered  in  order 
that  a  small  adjusting  wedge  could  be  placed  in  it  to  strain 
the  SMA  wire  to  a  specific  elongation  after  the  wire  was 
wound  around  the  small  teeth  of  the  L-shape  beam.   The  author 
used  compressed  air  to  blow  the  metal  crumbs  away  after 
trimming  and  sanding. 
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Figure  3.3.   Specimen  pre -curing  bagging  tool 
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The  author  used  a  0.02 "-thin  aluminum  separator  to 
separate  the  upper  and  lower  layers  of  SMA  wires,  as 
depicted  in  Figure  3.2  and  Figure  3.3. 

The  caul  plate  must  be  made  so  that  its  edge  could  just 
fit  against  the  L-shape  beam  to  prevent  cracking  the  tooth 
roots  of  the  L-shape  beam  by  excessive  wire  tension. 

Finally,   a  thin  sheet  of  mylar  was  placed  on  top  of  the 
whole  frame  and  a  release  agent  was  sprayed  on  the  mylar. 
Meanwhile,   the  whole  frame  and  the  SMA  wires  were  cleaned 
with  acetone.  Then,  a  thick  sheet  of  mylar  was  placed 
beneath  the  whole  assembly.  The  Allen  screw  shown  in  Figure 
3.2  can  be  adjusted  to  strain  the  wires  which  were  ready  for 
the  embedding  process. 

3 -l-S  First  Batch  Specimen  Fabrication 

The  prepreg  of  GF/Ep  uni -directional  tape  must  be 
cleaned   [9] .  Any  dirt,  dust,  or  other  foreign  objects  which 
come  in  contact  with  prepreg  tape  will  cause  voids  within 
the  hybrid  composite  laminates.  Therefore,   the  base  plate 
and  wire-straining  frame  should  be  as  clean  as  possible.  The 
tool  plate  and  the  caul  plate  had  to  be  sanded  with  number 
200  or  finer  sandpaper  until  the  surface  was  smooth  with  a 
dull  shine. 

The  author  covered  the  outside  inch  of  the  base  plate 
with  tissue  or  a  wood  frame  before  applying  the  mold  release 
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agent.  It  is  not  desirable  to  have  mold  release  agents  on 
the  outer  inch  of  the  base  plate,   for  this  area  will  be  used 
to  seal  the  vacuum  bag  on  the  base  plate.  The  mold  release 
agent  was  sprayed  onto  the  base  tool  plate,  both  sides  of 
the  caul  plate,   the  aluminum  dams,  and  wire -straining  frame. 
All  of  the  sprayed  pieces  were  placed  in  an  oven  and  baked 
at  250°F   (120°C)    for  approximately  two  hours. 

Initially,  before  straining  the  wires  between  the 
opposite  L-shapes,  wire  was  wound  on  a  hair-comb  set  (the 
same  tooth  size  as  the  L-shape  beams  in  Figure  3.2) .  Then, 
laid  up  the  first  two  layers  of  GF/EP  on  thin  mylar  and  the 
base  plate,  using  a  lamina  cutting  tool  specially  made  to 
cut  a  lamina  of  a  definite  size.  The  wound  wire  was  trans- 
ferred from  the  hair  comb  set  to  the  L-shape  beam  on  top  of 
the  first  two  layers.   Put  the  other  six  laminae  on  top  of 
the  wires.   Transferred  the  upper  layer  SMA  wires  from 
another  hair-comb  set  to  the  L-shape  beam  and  placed  the 
last  two  layers  of  laminae  in  position.  This  completed  the 
lay-ups  of  the  SMA  hybrid  composite  laminate. 

Next,   the  adhered  instruments,   such  as  the  embedded 
thermocouple  head  or  other  connecting  measuring  devices, 
were  cleaned  with  acetone  and  then  calibrated. 

Finally,   bagging  materials  were  wrapped  on  top  of  the 
wire-straining  frame  with  the  hybrid  laminates  for  curing 
(Figure  3.3). 
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The  fabricated  hybrid  laminate  is  a  20cm  x  15cm 
rectangular  laminate   [0/90/SMA/0/90] ^   (Figure  3.1),  which 
then  can  be  cut  to  the  assigned  specimen  size.  To  do  this: 

•The  SMA  wire  surface  must  be  sanded   (2  00  or  finer 
grade)   and  cleaned  with  acetone  after. 

•A  thermal -insulated  coat  may  be  required  in  order  to 
shorten  activation  time  for  the  SMA  wire  actuators. 

•Spray  release  agent  on  the  wires  outside  of  the 
laminate  for  easy  peeling  off  the  bagging  materials. 

•A  total  of  four  specimens  can  be  made.  See  Figure  3.1. 

3 • 1 -4 — Second  Batch  Specimen  Fabrication 

It  was  decided  that  the  wire  straining  frame  for  the 
first  batch  of  specimens  was  not  refined  enough.  Therefore, 
a  new  wire  stretching  tool  was  devised  for  more  efficient 
fabrications.  The  new  stretching  tool  and  the  new  base 
plate/dam  set  for  composite  laminates  are  shown  in  Figure 
3.4. 

3.2     SpecimPH  Testing  Plan 

(1)     The  following  items  needed  to  be  tested,   either  to 
verify  the  provided  data   [8]   or  to  get  more  test 
information  for  the  transient  thermal  response. 


•  Measuring  Nitinol  electric  resistivity  under  the 
temperature  ranging  from  60°F  to  22  0°F, 

•  Measuring  the  heat  conductivity  of  Nitinol  and  the 
whole  specimen   (which  required  a  calorimeter) 

The  laminate  beam  must  be  clamped  or  simply  supported 
at  two  ends.  The  modal  test  for  the  cantilever  beam 
will  be  carried  out  for  the  first  batch,   type  1  only. 
Because  the  free  end  will  be  connected  with  electric 
jumper  heads,   the  test  data  may  not  be  accurate  for  the 
cantilever  beam  . 
The  author's  test: 
•The  natural  frequency 

A.  SMA  hybrid  laminated  beam. 

B.  Gf/Ep  laminated  beam  (no  wires) . 

C.  SMA  hybrid  beam  with  electric  current  at 
steady  state  beam  surface  temperature  reaches 
austenite  phase. 

•Damping  coefficient  with  the  above  A,  B,   C  conditions; 

Perkins   [3]   shows  that  SMA  has  lower  electric 

resistance  at  -60°C. 
The  clamped  support  must  be  thermally  insulated  to 
assure  that  there  is  a  minimum  heat  loss  during  the 
heating  process.  The  insulation  material  may  be 
polyurethane  or  some  type  of  refractory  felt.  Another 
method  is  to  apply  an  insulated  pad  while  the  hybrid 
laminate  is  activated,   then  remove  the  insulated  pad  at 
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the  deactivation  process  or  use  airflow  to  cool  down 
the  hybrid  laminate  surface  so  the  SMA  wires  can  return 
to  the  martensite  phase. 
(4)     Electric  current  power  supply: 

The  power  supply  should  be  impulse  electric  current, 
machine  controlled  manually.  The  diagram  of  Figure  2.1 
shows  the  equipment  of  the  wire  connector,  circuitry, 
the  thermocouple,  and  the  FFT  dynamic  response  test 
devices   [4-7]  . 

 Facilities  for  Specimen  Fabrication.  Testing 

The  facilities 

(1)  Composite  parts  fabrication  autoclave. 

(2)  All  the  tooling  materials  and  f iberglass/epoxy  prepregs 
are  available  in  shops. 

(3)  SMA  wire   (300  feet  total  for  first  and  second  batch 
tests  of  SMA  hybrid  laminates) . 

(4)  Testing  the  natural  frequency,   and  the  damping  (the 
data  acquisition  system  in  the  structure  laboratory) . 

(5)  Testing  electric  resistivity  (a  common  multimeter) . 

(6)  An  impulse  electric  current  generator,   to  activate  SMA 
wires . 

(7)  Electric  wire  connecting  devices,   co-axial  cables. 
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Figure  3.4.     The  L-shape  beams  with  teeth,   the  base 
plate/Dam  set,   and  the  stretching  tool 
at  the  far  end. 


CHAPTER  4 

EXPERIMENT,   ANALYSIS  OF  SPECIMEN  TRANSIENT 
THERMAL  RESPONSE 

4^1 — Analysis  of  Transient  Thermal  Response 

The  simulated  thermal  model  is  described  in  Figure  2.2. 
The  X-axis  is  the  width,  and  the  y-axis  is  the  depth  of  the 
SMA  hybrid  beam.  Due  to  the  symmetric  wire  and  laminate 
combination,  the  small  cell  depicted  in  Figure  2.2  was  used 
to  verify  the  transient  thermal  response  given  by  Hawkins  et 
al.  [6]. 

The  simulated  thermal  model  is  used  in  predicting  the 
temperature  distribution  (transient  thermal  response)   in  the 
said  single  cell  of  specimen  type  1.  By  adjusting  the  DC 
power  applied  to  the  specimen,  author  could  compare  the 
ANSYS  analysis  results  with  the  experimental  results. 

Using  the  ANSYS  thermal  analysis  package,  the  author 
obtained  the  analysis  results  and  printouts  of   (case  la   [6] ) 
model  and  the  specimen  type  1  model.  Diagrams  of  these  plots 
of  temperature  distribution  and  the  temperature -time  can  be 
seen  in  Figures  4.1,   4.2,   and  4.3. 

From  the  analysis  results  and  plots,   it  was  verified 
that  the  temperature  inside  the  SMA  hybrid  composite 
laminate,  as  long  as  the  laminate  is  thin  compare  to  the 
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beam  length,  was  constant  and  uniform  at  each  time  step  (the 
temperature  differences  in  diagrams  were  not  much  off  from 
the  SMA  to  the  beam  surface) . 

4^l2  Experiment  of  Transient  Thermal  Response 

The  specimen  for  thermal  response  experiment  was  made 
from  fiberglass  prepreg  with  embedded  SMA  wires,  which  is 
described  in  the  previous  section.  The  facilities  of  testing 
were 

•  A  K  type   (TT-K-36  SLE  thermocouple  wire,   spot  welded 
under  1.3  watt  at  the  end),  a  SMP  male  connector (#1 
style  no. #10),   a  digital  thermometer   (model  HHM59--K 
type)\   and  a  thermocouple  socket. 

•  A  DC  power  source  and  a  Triple  output  DC  power  supply 
1660   (BK    Precision  Co.)   were  the  power  source  for  the 
activation  of  SMA  wires. 

The  analysis  tool  for  the  transient  thermal  response  is 
ANSYS.  The  geometry,  density,   and  the  thermal  constants  are 
listed  as  follows: 

p  =  404.3  lb„/ft^ 

(SMA  wire) 

1  =  0.005"   =  1.27  X  10-"  m  -  '        (SMA  wire)  [8] 

Pe=  82  X  10-«  ohm-cm  (SMA  wire)  [8] 

kx=  0.17  watt/cm/°c   (thermal  conductivity  of  SMA  wire) 

(   9.82  BTU/hr/ft/°F)  [5] 
Pc=  112.3  lb„/ft^  (density  of  composite)  [8] 

kcy=kcx=  k^/lS  (thermal  conductivity  of  composite)  [8] 
Cp=  1757  joul/kg/°k  =.3776  BTU/lb„/°F  (composite)  [5] 
h  =  29.71  watt/mV°k  =5.232  BTU/ftV°F   (composite)  [5] 


'Omega  Technology,   P.O.  Box  4047,   Stamford,   CT  06907-0047 
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The  experiment  on  SMA  wire  electric  resistance  to 
verify  the  electric  resistivity  provided  by  [8]   and  the 
connectivity  of  the  wires  inside  the  specimen  was  done  at 
the  laboratory  (A  common  multimeter  was  used  to  measure  the 
electric  resistance) .  The  specimen  was  the  Type  1  plate  with 
20  SMA  wires  each  on  the  upper  layer  and  lower  layers.  Both 
the  upper  layer  and  lower  layer  wires  at  the  first  entry  of 
the  specimen  were  crimped.  The  cross  section  of  the  wire  was 
2A  where  A  is  a  single  wire  cross-sectional  area. 

The  calculated  resistance 
Pe  =  76  X  10"^  ohm  cm         (martensite)  [8] 
A  =  TT  X   (0.005  X  2.54   )'  =  5.067  x  10""  cm^ 
R    =  Pe  /2A  =  0.075  ohm/cm 

There  were  20  sets  of  wires  in  the  specimen,  each  set 
is  22cm  long.  Therefore,   R,„,^i  =  0.075  x  20  x  22  =  33  ohm; 
the  specimen  test  result  is  34  ohm;  and  the  error  between 
the  two  data  is  1/34  or  2.9%.  The  result  showed  good 
agreement  with  calculated  electric  resistance. 

The  model  dimension  was  measured  from  the  specimen.  The 
first  batch  test  specimen   (TYPE  1,   SMA  wire  Vf=  3.27%),  and 
related  data  are  listed  in  Table  4.1.  The  photograph  of  a 
real  specimen  is  shown  in  Figure  4.5. 
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Table  4 . 1    Specimen  data 


1    TYPE  1 

1.13"x7.88"   ,  overstrained  %  SMA  wires 
[0/90/SMA/0/90/90/0/SMA/90/0]  FG/EP 

TYPE  1 
(base) 

1.13"x7.88"   ,   no  SMA  wire 
[0/90/0/90/90/0/90/0]  FG/EP 

The  input  data  to  generate  the  ANSYS  plot  in  Figure  4.1 
is  listed  in  page  39,   and  the  trasient  thermal  response 
test/analysis  comparison  is  in  the  Figure  '4.4. 
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Figure  4.1.  The  temperature  distribution  from  ANSYS  result 
for  (case  la   [6])   at  t=  0.02  hr 


A  =  a/1   =  5.0,    B  =  b/1  =  2.5,    1   =  0.005" 
A  =  width,  B  =  depth 
k^=  9.82  Btu/hr/ft/°F  (SMA) 
p  =  404  Ibm/ft^  (SMA) 
k^=  0.655  Btu/hr/ft/°F  (composite) 
p  =  112.3  Ibm/ft^ 
c  =   .3776  Btu/lbm/°F 
T  =  70°F 

h  =  5.232  Btu/hr/ftV°F  (air  convection) 
q  =  14471730  Btu/hr         (heat  generation) 
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When  blowing  cool  air  across  the  beam  surface,  and 
h  =5232  Btu/hr/ftV°F  (assumed) .  At  t=0.0005  hr=2  sec, 
temperature  drops  fast . 

Figure  4.2.  The  forced  cooling  of   (case  la   [6])   from  ANSYS 
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Figure  4.3.   Thermal  response  of  Type  1  specimen  (first 
batch),   t=  0.03hr  ,    I  =  0.185  amp. 
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Figure  4.4.     Thermal  response  of  specimen   (type  1),  first 
batch.   I  =  0.095  amp. 


The  first  batch  specimen  (type  1)  showing 
left  end  connected  with  wires  for  electric 
current . 


CHAPTER  5 

CONSTITUTIVE  EQUATIONS  OF  SMA  HYBRID 
COMPOSITE  LAMINATE 


S^J^ — Stress /Strain  in  SMA  Hybrid  Composites 

It  is  known  that  the  conventional  thin  plate  theory  has 
been  useful  for  predicting  the  stress/strain  relations  of 
the  composite  plate.  The  hygrothermal  effects  can  also  be 
calculated  by  using  the  conventional  theories  of  composite 
plates  to  obtain  ply  strains  and  stresses,   then,   employing  a 
failure  criterion  to  determine  whether  .or  not  failure  will 
occur  in  a  ply.  The  composite  structure  could  be  graphite 
/epoxy-Fiberglass/epoxy  or  combined  Kevlar-graphite/epoxy 
laminate  composed  of  orthotropic  plies  forming  a  homogeneous 
solid  plate  or  beam. 

The  SMA  hybrid  composite  is  a  totally  different 
structure  possessing  composite  plies  with  buried  SMA  wires. 
The  SMA  wire  diameter  is  almost  25  to  30  times  that  of  the ' 
reinforced  fiber  in  a  composite  ply.  By  using  the  SMA  wire, 
a  new  type  of  hybrid  structure  is  formed  with  highly 
inhomogeneous  properties  and  the  SMA  wires  are  often 
considered  as  inclusions  in  a  composite  structure.  The  SMA 
wires  tend  to  shrink  back  to  the  austenitR  ph;..c,^  ^hen 
activated  by  heating,  which  is  why  they  must  be  strained 
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before  curing.  Thus,   in  this  chapter,   the  stress/strain 
relation  is  still  based  on  known  theories,  but  the  recovery 
stress  will  be  emphasized  in  the  following  sections.  The 
overall  hybrid  composite  stiffness  changes  at  different 
activation  temperatures.  At  the  activation  temperature,  the 
recovery  stresses  in  SMA  wires  can  create  internal  stresses 
at  the  interface  of  SMA  and  composite  materials. 

5.1.1     Thermal  Strain  of  Laminate 

The  general  expression  of  the  strain  in  each  lamina  is 
given  by 


e  (t)  =  (  r(t)  -  r„  ) 

r„  =  initial  temperature  before  activation 
r(t)  =  transient  temperature   after  activation. 

According  to  the  experimental  results  presented  in 
Chapter  4,   the  temperature  inside  the  hybrid  laminate  is 
fairly  uniform  at  any  time  step.  The  strain  of  the  plate 
containing  the  initial  strain  and  curvature  terms  is 
represented  by  e,^ .  The  subscript  ij  signifies  the  element 
in  the  column  vector.  The  expression  E^^"  is  the  mechanical 
strain  of  the  Kth  layer  which  can  produce  stresses  in  each 


ply,  and  a^^  is  the  thermal  expansion  coefficient  of  the  K- 
th  ply,  which  varies  at  different  ply  orientations.  We  can 
write   (5,1)   in  matrix  form  as  in  (5.2)    (see  Agarwal  [10]). 


(5.2) 


"ex 

+ 

ttx 

Bit) 

^  y 

+ 

6  it) 

v" 

Ixy 

K 

+ 

^>^^_ 

dit) 

The  extreme  right  bracket  in  (5.2)   is  the  combined 
effect  of  recovery  stress  and  thermal  stress  in  a  single 
lamina  when  SMA  wires  are  embedded.  We  consider  the  recovery 
stress  and  induced  ply  stress  in  a  global  view  at  the 
moment.  Time  is  represented  by  (t) . 


=   [  Q]k 


Yxy 


I  Q  Vk  =   [  T  ]^  [  Q'  ]      f  T  ]  (5.3) 


[  T  ] =  Transformation  matrix 

^   ^     ^  K  ~    ^^Li  Et,  Glt,  .  .  . ) 

From  reference  [7] 


A  A 
"A        ^  A 


Eg  =SMA  Young' s  modulus 

B„  =  comp.  lamina  modulus 

flg  =  SMA  cross,  area 

\  =  comp.  lamina  cross .area 


(5.4) 
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~  a    F    f  f^^p     '  -^-"T  =  transverse  modulus  (5.5) 


vIt  =  +  v^:^  (5.6) 


Es:  The  SMA  Young's  modulus  should  be  time  dependent 
Es(t).   It  is  not  linear  with  respect  to  temperature  from 
raartensite  to  austenite  phase,  although  it  is  linear  from 
100°F  to  120°F  during  activation  [1,2]. 

The  above  equation  of  lamina  transverse  stiffness 
properties  could  be  modeled,  according  to  Halpin-Tsai, 
Agarwal    [10] . 

S.1.2  SMA  Hybrid  Laminate  Strain 

Therefore,   from   (5.2)   to   (5.7)   we  get 


K 

Ox 

el' 

k'' 

rh/2 
-h/2 

dz  =  [  A'  ] 

+  [  B-  ] 

<• 

(5.8) 

K 

47 


X 

W.. 


I- 


h/2 
h/2 


.M 


Z  dz  =   [  B*  ] 


y 


+  [  D-  ] 


(5.9) 


The  extreme  right-hand  term  of   (5.8)   contains  the 
thermal /recovery  in-plane  force  N/*  and  other  terms;  the 
extreme  right  hand  term  of   (5.9)   contains  the  thermal/ 
recovery  moment  M/*  and  other  terms.  The  M^^*,  My^*,   and  M^^* 
cease  to  exist  if  symmetric  SMA  lay-ups  are  formed. 


[A']  = 


rh/2  — 


[Q]  K  dz 

J  h/2 

[B']=\  [Q]\zdz 

J -h/2 

[Q]k  z^  dz 

J -h/2 


(5.10) 


ex 

+ 

+ 

A' 

B'' 

+ 

B' 

D' 

+ 

Mx 

My 

+ 

My 

+ 

M^ 

(5.11) 
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When  the  in-plane  force  and  moment   (external  load)  plus 
the  thermal /recovery  force  and  the  boundary  loads  (SMA 
activation)   as  suggested  by  Agarwal   [10]   is  applied  on 
laminate,  we  are  able  to  find  the  laminate  strain  as 
following . 


y 

y  xy 


A'  B' 
B*  D' 


-1 


+ 

+ 

+ 

+ 

Mr 

My 

+ 

+ 

ihybr.  laminate  strain) -^2) 


Note  that  in  the  column  vector  of  previous  equations  or 
in  the  later  equations,   the  superscript  *  of  each  term 
represents  laminate  properties  and  the  thermal /recovery 
stress  change  with  respect  to  temperature.  The  temperature 
changes  at  each  time  step  when  the  SMA  wire  in  the  laminate 
is  activated  or  deactivated. 


°y 


(t)    =   [  Q]'k 


«x' 

^y 

+  z[Q]k 

"y 

(5.13) 


Therefore,  due  to  the  laminate  strain,   the  Kth  layer 
stress  of   (5.13)   is  primarily  attributed  to  equation  (5.14), 
which  is  the  stress  resultant  of  a  1-D  beam.  The  symbols 
used  in  equation  (5.14)   are  similar  to  those  of  equation 
(10)   of  Rogers  et  al.  [7]. 


Nr  =  (I-XJ  £  a["'  v<"'(z„.,  -  z„) 
-  n-Kj  E  ai"'(l-v<"')  (z„.,  -  z„) 

N 


(5.14) 


The  N  layers  of  hybrid  laminate  have  recovery  stresses 
which  are  intermingled  with  thermal  stresses.  With  the 
additional  boundary  stresses,   the  stress  resultants  can  be 
estimated  as  shown  in  equation  (10)   of  Rogers  et  al .  [7]. 

ol"^  =recovery  stress 

=nitinol  boundary  condition 
a<l"'  ^'Thermal  stress 

=  thermal  boundary  condition 

=Boundary  stress 
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The  thermal  boundary  condition  K„  =  0.77  is  for 
"rested"  SMA  hybrid  composite  beam,  and  the  S^4A  boundary 
condition  is       =  0.71  from  the  test  results  of  Rogers  et 
al.  [7]. 

The  term  in  equation   (5.14)   is  the  induced  in- 

plane  force  in  x-direction  of  a  beam  structure.  The  N/*  and 
the  N^'^*  will  be  presented  in  2-D  SMA  hybrid  plate  struc- 
ture,  if  there  are  no  bending  and  twist  moments  induced  by 
nonsymmetric  activations.  The  equation   (5.14)  must  be 
calculated  prior  to  being  inserted  in  equations   (5.11)  and 

(5.12)  along  with  the  applied  external  loads  N^,  Ny  and 

to  obtain  the  hybrid  laminate  strain.  Then,  we  obtain  ply 
stresses  at  different  time  steps  as  described  in  equation 

(5.13)  . 

The  said  expressions,   however,   are  not  adequate  to 
determine  what  happens  at  the  interface  between  each  single 
SMA  wire  and  the  composite.  Therefore,  another  approach  is 
necessary  to  get  N/*  and  other  terms.  This  includes  looking 
into  the  internal  stress  distribution  caused  by  SMA  recovery 
stresses,  which  will  be  discussed  in  the  next  section. 

The  in-plane  force,   such  as  the  N/*,   is  the  stiffening 
force  in  a  beam  structure.  Because  of  the  in-plane  force, 
the  stored  elastic  energy  of  a  beam  structure  can  increase 
tremendously. 
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 Internal  Stresses  of  SMA  Hybrid  Rf?am 

The  SMA  wires  generated  large  recovery  stresses  at 
activation,  yet  the  internal  stresses  induced  by  the 
recovery  stresses  could  not  be  found  by  conventional 
methods.  Therefore,   the  sole  objective  of  this  section  was 
to  find  compressive  stresses  and  shear  stresses  attributable 
to  the  SMA  wire  activation. 

Nomenclature 

C  =  cross-sectional  width  of  composite  strip 
rs=  SMA  wire  cross-section  radius 

do=  SMA  wire  free  shrinkage  length  (ref .  to  strain) 
aR=  SMA  recovery  stress   (uniform  at  cross  section) 
0^=  compressive  stress  of  composite  strip  (0^=0^^) 
T  =  shear  stress  of  composite  strip  (t=t„) 
To=  shear  stress  at  SMA,   composite  interface 


5 • 2 . 1  Stress  Distribution  at  Interface 

During  activation,   the  SMA  wires  tend  to  shrink  inside 
the  laminated  composite,  but  the  composite  prevents  the  SMA 
from  recovering  its  original  length  (L-d^)  .  The  stresses  at 
each  cross  section  tend  to  balance  out  one  another  under 
cantilevered  boundary  condifinn     as  Figure  5.1  shows. 
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Cantilever  bean 

///mm///////////// 


45/-45 

berglass/  epoxy 

1  or  2  layers  of  prepreg 


Figure  5.1.   3-D  view  of  SMA  hybrid  composite  beam  at 
activation  temperature 


The  SMA  hybrid  composite  strip   (C  x  C  x  L)   in  Figure 
5.1  contains  one  or  two  layers  of  composite  laminae   (in  this 
case,   fiber  orientation  is  along  the  SMA  wire) .  Effective 
material  properties  are  assumed  homogenous  but  orthotropic 
in  z,   r,    e  directions.   The  SMA  wires  with  r^  =0.005"  are  all 
uniformly  allocated  in  the  composite  beam;  the  two  SMA  wires 


in  Figure  5.1  are  representative  of  still  other  wires  not 
shown  in  the  diagram. 

Because  the  net  in-plane  force  should  be  zero  at  any 
cross  section   (z=Zo)  under  cantilever  boundary  condition, 
the  in-plane  forces  must  balance  each  other  between  the  SMA 
wire  and  the  composite  strip.  The  SMA  is  stretched  by  the 
composite  strip  and  the  composite  strip  is  compressed  due  to 
the  shear  stress  exerted  by  the  SMA  wire  longitudinal 
interface   (see  Figure  5.1).    (It  should  be  noted  that  the  in- 
plane  force  may  not  be  zero  under  clamped -clamped  or  simply- 
supported  boundary  conditions.) 

It  is  possible  to  use  stress  equilibrium  equations 
(equations  of  motion)   to  find  stresses  in  this  region. 
According  to  Malvern   [11,   1.1.2],   the  theories  of  elasticity 
and  plasticity,  which  are  based  on  the  concept  of  continuous 
material,   lead  to  quantitative  predictions.  These  agree 
closely  with  experience  over  a  wide  range  of  conditions. 
Since  the  strip  is  continuous  and  is  assumed  to  have  no 
voids  or  flaws  in  the  composite  strip,  we  can  use  the 
equation  of  motion   [11,  p. 668     (II. 4. CIO),    (II.4.C11)]  to 
solve  the  problem.  Practically,   Poisson's  ratio  could  be 
ignored  in  beams . 

Based  on  the  above  assumptions,   the  significant 
stresses  in  the  strip  region  are  ct^,  and  t„,  which  can  be 
seen  in  Figure  5.2.    (stress  a,,  will  be  dubbed        and  stress 
T„  will  be  called  t.)    In  the  Figure  5.1  and  5.2,   we  pick 


one  composite  strip  containing  one  SMA  wire,  as  a  force 
equilibrium  model.   In  this  model,   the  interface  between  the 
SMA  and  the  strip  is  shear  stress  t^^  caused  by  the  SMA  wire 
shrinking  effect   (recovery  stress) . 


Figure  5.2.     Stress  diagram  in  strip  at  activation 
temperature . 
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The  body  forces  and  the  inertia  terms  due  to  acce- 
leration are  negligible.  By  neglecting  the  stress  terms 
o'rr/   (^68'   TrO'   'T zB  >   the  equations  become 


~  =     °  (5.15) 

0  z 


5^rz  a  a,. 


3  r  a  z 


=  0  (5.16) 


At  this  point,   there  are  only  two  non- trivial  equations 
left.   Equation   (5.15)   shows  the  following:   t^,  =  f  (r)  ,  the 
shear  stress,   induced  by  the  SMA  recovery  stress,  is 
constant  along  the  z-axis,   and  it  is  soly  a  function  of  r. 

The  second  equation  (5.16)   can  also  be  derived  by  the 
stress  equilibrium  diagram   (Figure  5.2) 


+  da^  /dz  dz 


T„  +  dr^^/dr  dr 
k--  dz  -^1 
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set  T„=T 


(  T  +  dr/ar  dr) (  r+  dr)   dd  dz     +     t  rdd  dz 
(        +  da^/dz  dz)   rde  dr    -       rdd  dr 


d  d    T  T 

  —          +    =     0  (5.17) 

0  z  d  r  r 


Since  T  =  f(r),  Equation  (5.17)   can  be  written 


as 


(z,r)  =  (  -i£i£l  -  £l£l 


)    2  +  C„ 


(5.18) 


The  shear  stress  distribution  flri  (r,,)   is  induced  by 
the  SMA  recovery  stress  a,.  Therefore,  by  using  the  force 
balance  diagram  in  Figure  5.3,  we  are  able  to  find  a 
relation  between       and  t„. 


Figure  5.3.     Recovery  stress,   shear  stress  diagram 


Agarwal    [10,    (p. 72, 73)]   cites  the  most  quoted  theory  of 
stress  transfer,   the  shear  lag  analysis.  The  SMA  wire 
diameter   (0.025  cm)   is  much  bigger  than  the  E-glass  or  the 
S-glass  fiber  diameter   (0.001  cm  =  lO/xm)   Agarwal    [10,  p.  17] 
If  the  composite  strip  is  treated  as  homogeneously 
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reinforced  polymer  matrix  surrounding  the  SMA  wire,   the  wire 
behaves  as  a  fiber  in  the  described  shear  lag  principle. 
Also,   Piggott  et  al .    [12,    (p. 110, 111  Fig  2,3)]   showed  the 
shear  stress  exerted  on  the  matrix  by  the  fiber  under  simply 
supported  ends  can  be  used  as  a  base  model  applied  in  the 
force  diagram  shown  in  Figure  5.3, 


Or  7rr%  +  27rr3  T  dz     =   (        +  da^  )  irr^ 
d  CTr  2  t 


d  z  r3 
2  T 


z 


r. 


(5.19) 


(5.20) 


The       is  known;  it  is  dependent  on  temperature  T  and 
the  initial  pseudoelastic  strain   [2,3]   before  curing.  So 

£.Ro:LJS:  (0<k<l)  ,  where  k  is  the  SMA  wire  slipping 

coefficient   (depending  on  the  matrix  )   and  a^,  is  the 
recovery  stress.    (For  conservative  prediction,  assume 
k  =  1.) 

We  assume  that  t  =  t,,  which  must  be  constant  along  the 
z  direction  in  Eq   (5.15),   is  also  constant  along  the  r 
direction.    (Refer  to  Piggott    [12]   with  modified 
assumptions.)   This  assumption  is  the  most  critical,  because 
we  know  that  usually  it  is  the  case  that  r  decreases  as  r 
increases  toward  a  free  surface.   If  this  is  so,   the  equation 
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(5.18)  becomes 


=  -  ^2  +  Co  =  ajr,z}  (5.21) 


The  magnitude  of        (5.21)   is  half  of       at  r=r,  ,  and 
Co  is  zero  for  the  case  of  a  cantilever  beam.  Thus,  we  use 
the  force  balance  to  equalize  the  SMA  recovery  force  in  the 
composite  strip   (see  Figure  5.4).   From   (5.20),    (5.21),  = 
-Oj,  r^  /2r    with  a'   =  ojoj,  ,   r'=  r/r^ 

a'   =  -l/2r'  (  1   s  r'   ^  s/r  ) 


At     crossection  oF  Z  =  Zq 


Figure  5.4.     Compressive  stress  diagram  in 
a  strip 
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Therefore,  from  the  force  equilibrium 

5 

(^R     rl    =  J      27rr  dr 


I 


(5.22) 


2iT  Tg  dr 


It  follows  that  from  equations   (5.20),  (5.21) 

^  =  ^  (5.23) 
For  different  laminae  lay-up  angle,  with  the  SMA  wire 
oriented  as  in  Figure  5.1,   the  o,  may  alter  at  various  ply 
constructions   (Co  not  zero)  .  The       can  depend  on  other 
factors  in  the  realm  of  microstresses ,   such  as  fiber 
buckling  strength,   interlaminar  shear  strength,   and  matrix 
cracking  strength  of  the  composite. 

The  shear  stress,   t,   is  actually  not  constant  along  the 
r  direction.  As  the  r  increases,    t  should  decrease.    (For  a 
simplified  solution  at  the  macromechanical  level,   the  Eq 
(5.23)   and  the  o'  diagram  in  Figure  5.4  is  sufficient.) 

For  a  more  realistic  distribution  of  a,  in  the  strip, 
we  can  approximate  a,  linearly  along  the  r  direction  as  well 
as  the  2  direction  (Figure  5.5).  The  diagram  of  a'   in  Figure 
5.4  actually  contradicts  the  force  equilibrium  in  Eq  (5.22) 
if  we  include  the  rest  of  the  curve  of  a'  beyond  r=2r,, 
which  is  not  zero  even  at  r  =  oo.  Therefore,  we  tend  to 
single  out  the  composite  stress  diagram  at  each  small  strip 
(C  X  C  X  L)   as  Figure  5.5  describes. 
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I 

X 

Qt      Z    =  Z 


Figure  5.5.     Modified  a,  distribution 

With 


(x-rjr,  (x-r  )r  ^  ^°  (r.srsx)  (5.24) 


Again,  using  force  balance  at  z=2,  .  we  can  find  x  , 

a,  27rrdr  (r^srsx)  (5.22) 


From  (5.24) 


^        IX  r^)r^         (x-r  )  r. 
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6  r^^^  =  3       +  3  r,x  -  2y}  -  2x       -  Iv^'^ 


X  =  ^   (  v/ii"  -  1  )        =    2.3723  (5.25) 


Therefore,  we  substitute  (5.25)  back  to  (5.24)  and  get 


0^-  (  14£1Z  -  .7287  4  )  r.z„ 


r,  ■  ^  '    '0  -^0  (5.26) 


This  is  the  compressive  stress  at  z=Zo.  Putting  (5.26) 
back  into   (5.18)   to  redistribute  r  =  f (r) .  we  obtain 


(    1-7287  _    .7287  ^  j   _  _   d  f(r)  £(r) 

d  r      ~  (5.27) 


The   (5.27)   should  hold  at  any  z  =  z,  section.   It  is  a 
typical  O.D.E,  which  we  can  solve  for 
faCr)    :  homogeneous  solution 

f„(r)   =  c,/r 
fp(r)    :  particular  solution 

fp(r)   =   (rV  (4.1169  r,^)     -   .8284  r/r3)  r, 
with    B.C.    :         T     =  f(rj   =  To 
we  obtain 


■•'oSKtvs;' 
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7=  (1.6215(^«)-. 8644(^1)+. 2429  (J:)2)t„       (r.srsx)  (5.28) 


Checking  Eq   (5.28)   at     r=  2.3  723        (the  value  of  x)  ,  we 
have  T  =  0.0001,    j,  =  0,   and  r  =  at  .(checked) 

Finally, the  last  part  is  how  do  we  find  constant  Tq. 
From  Eq   (5.20)   and  Duerig   [2,  p. 220  and  p. 222,   Fig.   3]  for 
Ni-Ti  SMA,   the  initial  strain  is  6_i;  so        is  about  400  Mpa 
(58  Ksi) .  We  have  found  that  the  value  varies  with  different 
SMA  materials,   so,   from  (5.20), 

58000  =  2   Tjr,     L  (r,  =  0.005") 

To     =  145  /  L    psi.  (5.29) 
Mohamed  [13]  also  showed  a  formulation  for  calculating  the 
SMA  recovery  stress  as  follows 

T  -  A 

o^{T,e)  =  ol'P  (  1  -  exp'"'^'  )    (  1  -  exp-"'  )  (5.30) 

"^^^^  (7r  =  recovery  stress  of  SMA 

a^tp^=  yield  stress  of  SMA  at  Af 

K  =  constant  analog  to  chemical  reaction  in 

austenite  phase  transform. 
Ag  =  austenite  starting  temperature 
Af  =  austenite  final  temperature 
N    =  number  of  crystallographically 
equivalent  martensite  variants. 
£    =  initial  given  strain 
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The  features  of  the  recovery  stress  of  (5.30)  are: 

(i)  For  a  given  strain,   £,        is  sensitive  to  temperature 
only  within  the  range        -  Af  . 

(ii)  If  (Tr  achieves  its  maximum  value  at  the  Af  temperature, 
the  equation  (5.30)  becomes 

max.  =  a^/^  (  i  -  exp"'''  )  (5.31) 

(iii)  For  temperatures  higher  than  Af,        becomes  sensitive 
only  to  initial  strain  e . 

The  above  can  refer  to  Duerig  et  al .    [2]   in  which  the  value 
of  Ne  can  be  found  through  cross-checking  the  above 
reference   [2,   p. 220  -  p. 222,   and  Perkins,  3]. 

^■'2-.2 — Micromechanics  and  Internal  Stre.g.qpg 

If  we  combine   (5.29)  with  (5.26)   and  (5.28),  we  get  a, 
and  T  in  the  hybrid  composite.  Then,  we  add  them  to  the 
laminate  hygrothermal  residual  stresses  and  in-plane 
stresses  due  to  the  static  and  dynamic  loading.   In  this 
moment,   the  ply  stresses  derived  from  classic  laminate 
theory  can  be  calculated. 

Finally,   from  the  work  of  Chamis    [14]   and  Guynn   [15] , 
microstress  level  can  be  checked  at  each  layer  of  the  hybrid 
laminate  to  see  if  any  of  them  exceeds  the  allowable  fiber 
strength  and  the  allowable  matrix  strength.   The  formulations 
in  Chamis   [14]   can  predict  the  microstresses  of  tension, 
compression,   and  shear  near  the  fiber/matrix  or  the 
matrix/matrix  interface,. 
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The  highest  a^,        at   (z  =  L  .  r=r„)   occurred  in  the 
cantilever  hybrid  beam.  The  compressive  stress  in  the 
composite  strips  may  be  redistributed  if  oJr^,h)  exceeds 
the  compressive  strength  ol^^  of  the  composite   (Agarwal  [10, 
p. 49-52,  Eq  (2.50)).   In  this  case,   the  clamped  boundary  may 
provide  more  constraint  in  redistribution,  and  the 
microstresses  at  this  region  can  be  critical  during  the  SMA 
hybrid  composite  beam  activation.  Recovery  stress  plays  an 
important  role  in  this  specific  micro-mechanical  analysis 
(Figure  5.6).  Different  boundary  conditions  for  the  hybrid 
composite  beam  have  different       distributions  along  the  z 
axis   (as  shown  in  Figure  5.7). 


repeated 

0°  edge  of  beam)   (io  edge  of  beam) 


—  C/2 


At  Z  =  L 


Figure  5.6.     Clamped  boundary  diagram 
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The  Full  clamp- clamp  condition   (Figure  5.7)   creates  the 
term  Njl  (section  5.1.2,  Equation  (5.18))  which  is  the  in- 
plane  force  at  a  particular  cross  section  of  a  clamped- 
clamped  beam. 


Cantelevered 


Simply  support 


Full  clamp-clamp 
(create  inplane  force) 


Part  clamp-clamp 
(create  Inplane  force) 


L/2      L  Z 


(Jr. 


Z  =  Z, 


CT. 


z  =  z, 


1 

Oc  r 

z  =  z, 


5 


z  =  z, 


Figure  5.7.     Four  examples  of  internal  stres 
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The  full  clamp-clamp  condition  in  Figure  5.7  is       =  Cq 
or       =  0  ,   T  =  To/r    or  t  =  0;  and  the  Cq  of  Eq  (5.21) 
depends  on  the  boundary  clamping  device  and  the  thermal 
conditions . 

 Dynamic  Model  of  SMA  Hybrid  Laminate 

The  concept  of  using  SMA  in  reinforced  laminates  was 
found  in  Rogers   [1],  Saunders  et  al .    [5],  and  Rogers  et  al . 
[7] .  The  benefits  of  stiffening  part  of  a  structure  are 
obvious  in  critical  situations  where  resonances  or  acoustic 
noises  impinge  on  the  structure.  SMA  has  the  unique 
properties  of  promoting  modulus  variation  and  change 
recovery  stresses. 

5.3.1     General  Dynamic  Model 

The  governing  equation  of  a  symmetric  SMA  hybrid 
laminated  plate  can  be  derived  by  using  the  Ritz  method  or 
the  energy  variational  principle,  both  mentioned  in  Saunders 
et  al.    [5]   and  Rogers  et  al .    [7].  we  assume  that  the 
structure  is  excited  by  a  harmonic  load,  the  general 
governing  equation  is  described  as  : 

dxdY'      ''ay^      ^aF  ^"ay^ 
+  2N'         +  r      ^  .  ,      .    ■  . 
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where 

W  =  deflection  function  W{x,y) 

=  damping  coeficient 
p  =  mass  per  unit  area 

2(x,y)  e-'"'^  =  external  forcing  function 
K  =        =  SMA  recovery  stress  resultant 

This  flexural  equation  can  be  expanded  to  handle  the 
response  of  an  asymmetric  laminated  plate.  The  W{x,y)   can  be 
solved  either  by  employing  FEM  or  the  Galerkin  method,  the 
choice  depending  on  the  boundary  conditions. 

5-3.2  Dynamic  Model -Qrthotropic  Laminate  (APT^ 

The  orthotropic  laminates  are  commonly  used  in  many 
composite  structures,   and  its  SMA  hybrid  laminate  can 
perform  as  active  properties  tuning   (APT)   elements.  The 
hybrid  structures  are  usually  cantilever  plates  with  linear 
vibration  characteristics.  Their  flexural  stiffness  D,^  is 
the  only  temperature  dependent  term  at  the  SMA  activation 
and  there  are  no  in-plane  forces. 

dX^  dY^  dY^  dt^ 


^•^•3  Dynamic  Model -Qrthotropif^  Laminate  (ASET) 

The  governing  equation  of  this  type  of  structure  is 
always  under  specially  prescribed  boundary  condition  such 
hinged-simply  support,   clamped-clamped,   and  hinged- clamped 
conditions . 
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This  means  that 


For  example,   in  a  1-D    hybrid  composite  beam  with  a  hinged 
simply- supported  condition,   the  recovery  stress  resultant  of 
SMA  is  not  zero  and  the  equation  of  motion  can  be  described 


"  ax*       dx^     ^  d  t2 


The  recovery  stress  resultant,   combined  effects  of  the 
thermal  and  prestressed  condition,   is  present  in  the 
equation  while  the  SMA  is  activated,  and  we  assume 

W  =  W{x)  e^"'^ 

Then  the  natural  frequency  of  the  mth  mode  is 


L 


Apparently,  the  natural  frequency  w„  is  higher  than 
that  of  a  beam  which  is  not  activated.  The  first  or  higher 
modes  of  natural  frequencies  can  be  shifted  by  heating  the 
SMA  wires  in  the  hybrid  composite  beam  to  avoid  resonance. 
It  is  expected  that  the  vibration  can  be  subdued  by  the 
activation  or  deactivation  in  SMA  hybrid  laminates. 


CHAPTER  6 
EXPERIMENTS  OF  SMA  HYBRID  LAMINATES 
MODAL  RESPONSE 


6.1     Experiment  Facilities 
The  equipment  used  in  the  testing  are  : 

•  LeCroy  8013A  waveform  recorder,  a  data  acquisition  device 
for  dynamic  response  of  structures,  LeCroy  Research 
System  Corp. 

•  Hammer  with  force  transducer,  an  impulse  load  on  beams. 

•  Eddy  current  probe  to  measure  beam  vibration  responses. 

•  K  type(TT-K-36  SLE)  thermocouple  wire  spot  welded  to  the 
end  to  an  SMP  male  connector  (#1  style  no.  10) .  Digital 
thermometers  Model  HHM59,  K  type  thermocouple  socket. 

•  IBM  PC-286  compatible,  equipped  with  CATALYST/AS YST 
software  which  analyze  the  signal  transmitted  from  the 
LeCroy  8013A  wave  form  recorder.  It  creates  dynamic 
response  data  in  the  time  and  frequency  domains. 

•  DC  power  source.  Triple  Output  DC  Power  Supply  1660  (BK 
Precision  Co.)  provides  the  SMA  activation  power. 

»  A  16.5  cm-thick  granite  slab  as  a  foundation  on  which  the 
impulse  load  hammer,    eddy  current  probe,    and  all  test 
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OMEGA  Technologies  Company,  Box  4047,  Stamford,  CT  06907- 
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facilities  are  set.  There  are  several  vibration 
isolators,  made  of  hard  rubber,  placed  under  the  slab  to 
eliminate  external  noise  signals. 


Figure  6.1.     Experimental  facilities  includinq  the 
LeCroy  8013a. 


6.2     1st  Batch  of  Specimen  Test 
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The  purpose  of  the  test  is  to  verify  the  theoretical 
prediction  and  to  find  the  natural  frequency  and  damping 
factor  of  the  SMA  hybrid  laminate  beam.  The  first  batch  of 
test  specimens  has  two  types  of  SMA  hybrid  laminated  beams 
(Table  6.1),     as  well  as  two  base  laminated  beams  without 
SMA  wires  embedded.  We  will  compare  both  the  natural 
frequency  and  the  damping  factor  of  each  SMA  hybrid 
laminated  beam  to  its  base  laminated  beam  at  the  room 
temperature   (no  activation)   and  the  activation  temperature. 


Table  6.1     First  batch   (TYPE  1  ,  Vf=  3.27%  SMA 


j        TYPE  1 

1.13"x7.88",   overstrained  SMA  wires  1 
[0/90/SMA/0/90/90/0/SMA/90/0]  FG/EP 

TYPE  1 
(base) 

1.13"x7.88",   no  SMA  wire 

[0/90/0/90/90/0/90/0]  FG/EP 

(TYPE  2   ,   Vf=  1.63%  SMA  wire) 

TYPE  2 

l.l"x7.88"   ,  overstrained  SMA  wire 
[0/90/SMA/0/90/90/0/90/0]  FG/EP 

1  TYPE  2  (base) 

l.l"x7.88"    ,   no  SMA  wire 

[0/90/0/90/90/0/90/0]  FG/EP 

With  the  Type  1,  and  Type 
the  test  of  impulse  load,  only 
frequency  were  recorded  (Table 
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1  base  cantilever  beam  under 
the  first  mode  natural 
6.2)  . 


Table  6.2    Natural  Frequencies  of  cantilever  beams 


1*-           5.88"         ^  1 

TYPE  1 

TYPE  1  (base) 

I 

V 

Temp. 

1st 
mode 

Temp. 

1st  mode 

0. 

0  . 

R(75°F) 

42  HZ 

R(75°F) 

41.5  HZ 

.37 

amp 

24 

volt 

201°F 
(2  min) 

4  7  HZ 

/////// 

///////// 

The  difference  of  first  mode  natural  frequency  for  Type 
1  and  Type  1  base  is  not  much.  The  reason  may  lie  in  the 
canceling  effect  of  the  stiffened  wire  embedded  in  the 
laminate  minus  the  extra  mass  at  the  free  end  of  the  beam 
(crimped  metal  pieces  and  silver  paint  for  SMA  wire  electric 
connection)  under  cantilever  boundary  condition. 


The  natural  frequency  for  activated  Type  1  is  only 
slightly  higher  than  non-activated  Type  1  in  Table  6.2,  It 
appears  that  the  beam  has  a  very  low  SMA  volume  fraction 
(3 .27  %)  . 

The  results  of  Type  1  beams  under  impulse  load  and  Its 
first  mode  natural  frequencies  recorded  at  76.6°F  and  75°F 
as  well  as  the  activated  temperatures  at  200. °F,  313 °F  can 
be  seen  in  Tables  6.3  and  6.4. 


Table  6.3    Natural  frequency  of  clamped -clamped  beam 


1  cm    -»|      \*-  F      thermocouple  wire  (AL  foil) 

Vise  1    E  ©---E    vise  2  (17/12) 

|«-        14.5  cm  ^[ 

TYPE  1 

I 

V 

Temp. 

1st  mode 

0 

0 

76 . 6°F 

319.4  HZ 

0 

0 

76 . 6°F 

322.6  HZ 

*-  higher  clamp 
force 

200 . °F 

322.0  HZ 

Table  6.4    Natural  frequency  of  clamped-clamped  beam 


1  cm  -» 

vise  1 


|«-  F        thermocouple  wire  (AL  foil) 

vise  2  (18/12) 


14 . 5  cm 


(high  clamp  force) 


TYPE  1 

TYPE  Kbase) 

I 

V 

Temp. 

1st 
mode 

Temp 

1st 
mode 

0 

0 

75  .  °F 

336  HZ 

.  0095 

79°F 

307  HZ 

.  0186 

.52 

29v 

313°F 

324  HZ 

.  0088 

263° 
F 

300  HZ 

.  0323 

^   :  damping  factor 


hot  air  blowing 


The  conclusions  to  be  drawn  from  the  clamped-clamped 
beams   (Tables  6.3  and  6.4): 


Higher  clamping  force  gives  higher  natural  frequency. 
Activated  hybrid  laminate  beam  has  a  slightly  lower 
natural  frequency.  The  most  likely  reason  is  that  the 
SMA  wires  were  overstrained   (did  not  shrink  back  to  the 
original  austenite  phase  before  straining) ,   so  their 
recovery  stresses  could  not  overcome  the  thermally 
induced  beam  compression.  Therefore,  the  overall  in- 
plane  force  N/*  is  compressive.   Proper  choice  of  ply 
orientation  in  order  to  reduce  effective  thermal 
expansion  may  decrease  the  thermal  stresses. 
The  activated  hybrid  laminate  beam  has  a  lower  damping 
factor  than  its  non-activated  phase,  but  the  base  beam 
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(without  SMA  wires)   has  a  higher  damping  ratio  at 
higher  temperature  than  that  of  the  low  (room) 
temperature.  [16] 

It  should  be  noted  that  he  original  test  bed  in  the  lab 
was  designd  only  for  cantilever  beam  tests.  Some  insert 
threads  in  the  thick  aluminum  base  plate  would  need  to  be 
retapped  to  hold  the  second  vise  in  order  to  clamp  the  other 
end  of  the  beam  for  clamped- clamped  test.  The  insulation 
material  between  the  beam  and  the  vises  is  graphite  (carbon) 
felt  to  prevent  heat  conduction  from  beam  to  vises.  (The 
first  batch  test  specimen  did  not  have  graphite  felt 
insulator  at  the  clamping  location.) 

The  first  batch  specimens  were  not  satisfactory  and 
could  not  verify  that  an  activated  SMA  hybrid  laminated  beam 
had  a  higher  natural  frequency  than  a  non-activated  SMA 
hybrid  laminate  beam.  Some  specific  effects  were  found  in 
the  clamped- clamped  test   (Table  6.3)   that  showed  that  an 
activated  SMA  hybrid  laminated  beam  has  a  little  lower 
damping  factor  (^)  when  temperature  increases.  Compared  with 
the  base  laminated  beam  (no  SMA  wire) ,   the  damping  factor 
increases  as  temperature  increases   [16] .    (The  higher 
temperature  on  the  base  laminated  beam  is  imposed  by  forced 
hot  air.) 

The  first  mode  frequency  response  plots  of  the  data  from 
the  first  batch  specimens   (Table  6.4)   are  in  Figures  6.4, 
6.5,    6.6,   and  6.7. 
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6.3     Calculated  1st  Mode  Natural  Frequency 

Table  6.2  lists  the  experimental  data,  that  can  be  used 
to  verify  the  theory. 
(1)     Cantilever  beam:  Type  1  (base) 

First  approach  (simple  lumped  mass  approach) : 

1  =  5.88"   ,     t  =  0.076"    (8  layer  GF/EP) 
[9]     m  =  0.056  X  1.8/1.61  x  vol.   =  0.0354 
*  EI=  149.05  Ibf.in^ 

S  mm  a  =[  1^  /(  3  EI  )   ]   X  P 

so,   K  =    (  3  EI   )/  1^ 
03  =/(K/m)    =/(149.05  x  32.2  xl2)/(.0354  x  5.88^) 
=  374.47  rad/sec    =     59.6  HZ 


Second  approach  (closed  form  solution) : 
1  =  5.88",  t=0.076" 

p  =  0.056  X  1.8/1.61  X  t=0.005  Ib^/in^ 
*     EI=  149.05  Ibf.in^ 
[17]     CO  =  3.515  /(EI)    /  (pi") 

=  345.04  rad/sec  =     54.91  HZ 


Third  approach  (ANSYS) :  u  =  59.1  HZ 


*     GENC0M2:  a  self -developed  FORTRAN  code  for  general 
laminate  properties  calculations   (for  hybrid  laminate  also) 
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Comparing  the  second  approach  with  the  test  result (  w  = 
41.5  HZ),  we  see  that  the  error  is  about  24%.   It  clearly 
indicates  that  neither  the  clamping  device  nor  the 
fabricated  specimen  had  very  good  qualities,  or  there  are 
measured  errors.  Therefore,   further  improvements  in  the 
specimen  fabrication  and  testing  method  are  crucial  for 
future  study  on  sandwich  structures.  The  specimens  are  shown 
in  Figure  6.2. 

6.4     2nd  Batch  Specimen  Test 

The  lay-up  of  Fiberglass/epoxy  prepreg  with  the  SMA 
wires  is   [90/SMA/0/0/SMA/90] ;  the  0  degree  Fiberglass  fiber 
direction  is  aligned  to  that  of  the  SMA  wire  direction.  Just 
as  in  the  first  batch  specimens,   the  first  layer  is  20 
wires,  20  wires  down,  width  is  1.12",  and  volume  fraction  of 
the  SMA  in  laminate  is  6 .37% . 

The  specimen  is  shown  in  Figure  6.3.  The  test  results 
are  seen  in  Figure  6.8,   6.9,   and  6.10.  We  found  that  at  the 
activation  temperature  the  natural  frequencies  have  changed 
for  all  three  modes. 
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Figure  6 . 3 


Second  batch  specimen.  Type  1  (left) 


Response  of  SMA  hybrid  beam  (75°  F)  clamped-clamped 
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Figure  6.4.  First  mode  natural  frequency  of 
specimen 


Response  of  SMA  hybrid  beam  (313°F),  clamped-clamped 
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400  ■ 
350  - 


^     300  - 


0  200  400 


Frequency  (HZ) 


Figure  6.5.   First  mode  natural  frequency  of  Type  1 
(activated) 
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83 


Response  of  base  beam  (79°  F)  -  no  SMA  wire 


200  400 


Frequency  (HZ) 


Figure  6.6.  First  mode  natural  frequency  of  Type  1  (base) 
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Figure  6.7.   First  mode  natural  frequency  of  Type  1  (base) 
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Second  batch,  Type  1  (73°  F),  clamped-clamped 
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Figure  6.8.     Frequency  response  at  non-activation 


Second  batch,  Type  1  (140.6°  F),  clamped-clamped 


gure  6.9.     Frequency  response  at  activation 
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Figure  6.10.       Natural     frequency- -Temperature     plot  of 

second  batch  specimen 


CHAPTER  7 
VIBRATION  CONTROL  OF  SANDWICH  BEAMS 

Two  kinds  of  vibration  control  techniques  which  are 
currently  applied  to  reduce  the  problem.  One  is  passive 
vibration  control  technique,   the  other  is  active  vibration 
control  technique.  The  emerging  smart  materials  have  created 
other  ramifications  in  vibration  control  studies.  Chapter  1 
showed  the  unique  abilities  of  smart  materials  to  change 
material  bulk  modulus,   stiffness,  and  so  forth,  when 
activated  by  electric  current,  by  magnetic  field,   or  by 
heating.  The  ability  of  SMA  wire  to  change  stiffness  in  a 
SMA  hybrid  laminated  beam  has  been  described  in  Chapter  6. 
There,   the  test  results  showed  the  potential  of  the  hybrid 
material  to  shift  the  structure's  natural  frequencies.  Thus, 
the  resonant  condition  of  vibration  induced  by  impinging 
loads,  sinusoidal  loads,  etc.,  could  predictably  be  subdued 
by  activating  the  SMA  wires  embedded  in  the  composite 
structures.  However,   the  test  results  in  Chapter  6  also 
showed  lower  damping  coefficients  when  the  SMA  hybrid 
laminated  beam  was  fully  activated.  All  the  tests  of  SMA 
hybrid  laminated  beams  were  done  under  the  clamp-clamp 
boundary  condition.  A  simply- supported  beam  or  a  cantilever 
beam  is  often  chosen  as  a  candidate  in  practical  vibration 
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problems.  Still,   the  tests  showed  strong  evidence  of  a 
structure-stiffening  effect  and  the  changing  of  damping 
properties. 

In  practical  applications,   simply  supported  beams  and 
cantilever  beams  are  commonly  used.  When  a  beam  is  fixed  at 
two  ends  for  a  simply- supported  case,  or  a  helicopter  rotor 
blade  and  a  fan  blade  of  a  jet  engine  for  cantilever  beams, 
to  name  just  a  few.  The  simply- supported  and  cantilever 
beams  require  strength  and  bending  rigidity  to  serve  as 
sound  structures.  One  of  the  best  choices  for  a  simply- 
supported  beam  construction  is  a  sandwich  beam  .  There  are 
many  advantages  of  using  sandwich  structures:    (1)  sandwich 
structures  have  great  strength/weight  ratio,    (2)   they  have 
superior  flexral  rigidity/weight  ratio,    (3)   the  tailored 
composite  skins  perform  excellently  under  various  design 
requirements  and,    (4)   the  sandwich  core,  which  is  softer  and 
weighs  less,  can  be  made  out  of  good  damping  materials  in 
order  to  reduce  acoustic  loading,  or  other  applications.  The 
bonding  material  can  be  film  adhesive  and  spray-on  primer 
which  serve  the  purpose  of  join  sandwich  faces  with  core 
materials.  Thicker  core  can  be  cut  into  smaller  blocks  then 
strung  together  to  make  curved  sandwich  panels. 

The  candidate  sandwich  panel  in  this  chapter  is  a 
^^^^^"1   ^^^"^  whose  faces  are  E-glass/epoxy  containing  SMA 
"^^^^^  with  wire  volume  fraction  6.7%   (each  wi  rP  has  O.OT' 
diameter) .  This  adaptive  structure  may  be  activated  by 
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electric  current  or  by  a  direct  heating  process  while 
vibration  occurs. 

7.1    Vibration  Suppression  in  Structures 

Vibration  suppression  is  a  very  broad  subject  in  which 
many  techniques  are  available  to  achieve  this  purpose. 
Passive  vibration  control  is  a  technology  that  uses  extra 
mass-spring  attached  to  the  main  structures  or  viscoelastic 
damping  material  on  the  structures.   It  requires 
predetermined  dynamic  system  analysis  and  is  an  open  loop 
control  system.  Active  vibration  control  is  essentially 
closed  loop  control  to  suppress  mechanical  vibration.  It 
requires  the  effort  in  designing  feed  back  control  logic  and 
the  electronic  devices  in  structure  system  such  as  a  large 
space  station.  From  surveying  the  research  in  these  two 
control  systems,   the  author  found  that  the  SMA  hybrid 
sandwich  structure  possesses  both  abilities  in  passive 
vibration  controrl  and  active  vibration  control. 

7.1.1     Passive  Vibration  Control 

To  subdue  vibration,  one  could  use  an  added  mass  on  a 
vibrating  beam  as  a  typical  passive  control  scheme.   It  is  an 
undamped  dynamic  vibration  absorber  which,   according  to 
Thomson   [17]   and  Beards   [18] ,   can  reduce  structure  vibration 
but  extra  viscous  damping  will  have  to  be  added  to  reduce 
vibration  in  general  cases. 
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A  case  studied  by  Pierce  et  al .    [19]   uses  the  passive 
damping  principle  on  a  helicopter  rotor,  has  a  mass  M  and  a 
single  flap  pendulum  with  m  installed  on  the  rotor  to  absorb 
vibration.  This  set-up  would  cut  down  shear  force  transfer 
to  the  rotor  hub,  and  reduce  the  rotor  deflection  in  flight. 
The  mathematical  model  for  the  case  is  tedious  to  develop, 
and  the  method  required  trials  and  errors  to  locate  the  flap 
pendulum  on  the  rotor  to  achieve  ideal  results.   In  addition, 
in  real  structures,  the  design  for  the  installation  of  a 
flap  pendulum  on  a  helicopter  rotor  is  difficult.  An  extra 
viscous  damper  may  also  be  necessary  for  various  loading 
conditions,  not  to  mention  the  maintenance  for  the  rotor 
with  added  pendulum  in  the  long  run. 

In  viscoelastic  damping,   composite  materials  have 
better  damping  quality  than  most  monolithic  materials   [20] , 
[21] .  Short  fiber  composites  have  higher  damping  property 
than  unidirectional  fiber  composites   [22] ,    [23]   but  inferior 
stiffness  for  the  short  fiber  composite  is  inevitable.  They 
are  all  fiber  reinforced  polymer  materials,   and  the  polymer 
is  viscoelastic  material  which  provides  most  of  the  damping 
factor.  Therefore,  using  viscoelastic  damping  tape  is  more 
effective  for  increasing  the  damping  property  of  modern 
structures.  The  damping  tapes  in  structures  can  serve  as 
dampers  so  that  elastic  vibration  energy  can  be  dissipated 
as  heat    (Jones   [24] ,   Sankar   [25] ) .  Also,   in  a  large  space 
truss  whose  bar  elements  were  treated  with  segmented 
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constrained  damping  layers,  the  loss  factor  was  tremendously 
increased  --an  example  of  direct  evidence  of  elastic  energy 
dissipation  in  a  structure. 

The  added  mass  damper  could  increase  weight  and  not 
easy  to  apply  to  real  structures.  The  viscoelastic  damping 
treatment  is  easier  to  be  bonded  on  structure,  but  may  be 
peeled  off  by  external  disturbances.  The  sandwich  structure 
has  good  damping  core  to  be  used  for  passive  control  in 
noise  reduction  or  other  purposes.   If  SMA  hybrid  laminates 
can  combine  the  sandwich  core  become  a  property  tunable 
structure,  would  it  be  more  efficient  to  have  this  dual 
vibration  control  materials  in  one  structure? 

7.1.2     Active  Vibration  Control 

In  recent  years,   the  active  vibration  control  technique 
has  appeard  to  provide  more  options  to  control  mechanical 
vibrations.  An  example  is  the  alteration  of  the  stiffness  of 
several  bars  of  a  cantilever  truss  under  various  modes  of 
vibration  (Onoda   [26] ) .  Onoda  indicated  that  one  of  the 
control  laws  is  to  reduce  the  stiffness  or  increase  the 
stiffness  of  a  particular  member  of  the  truss  by  activating 
the  attached  piezoelectric  actuator  when  the  deflection  of 
that  member  exceeds  a  threshold  deformation.  So,   the  member 
operates  like  a  softener  to  absorb  vibrational  energy  or  a 
stiff ner  to  reduce  vibration.   This  unique  device  requires 
very  special  logic  control  circuitry,   and  the  actuator  tube 
connecting  the  adjacent  truss  members  also  demands  accurate 
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tolerance  control.  The  method  is  useful  in  space  structures 
and  many  mechanical  components. 

In  contrast  to  the  compicated  feed  back  control  system, 
the  SMA  hybrid  structure  could  be  modeled  as  a  simpler  feed 
back  control  element  with  the  passive  damping  core  material 
added  into  it . 

7.2     Equation  of  Motion  of  Sandwich  Panel 

In  this  section,   the  primary  goal  is  to  establish  the 
principles  of  mechanics  for  sandwich  panels  with  faces  made 
of  SMA  hybrid  composite  laminate.  A  large  flexible  structure 
like  a  wind  turbine,  at  random  wind  directions  plus  terrian 
variations,  may  face  low-frequency  vibrational  problems.  We 
are  assuming  that  this  flexible  wind  turbine  is  made  of  a 
sandwich  panel  which,   in  turn,  vibrates  with  large 
amplitude.  The  structure  is  under  elastic  deformation  when 
vibrates,  and  can  cause  large  amplitude  which  creates 
geometric  nonlinearity .  There  are  normal  or  in-plane 
stresses  in  the  core  with  constant  thickness.  The  only 
significant  stress  in  the  core  is  the  transverse  shear.  The 
face  skins  are  thin  compared  to  the  core  thickness  and  the 
upper  skin  is  identical  to  the  lower  skin,   and  the  face 
skins  can  only  take  the  in-plane  forces.  The  panel  is 
assumed  to  have  perfect  bond  between  the  skins  and  the  core. 
This  gives  the  most  general  case  of  kinematic  assumptions. 
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7-2. 1  Derivation  of  Equation  of  Motion 

Assume  this  general  case  of  a  cantilevered  sandwich 
panel,  only  one  face  is  activated,  erects  the  nonsymetric 
activation  model.  Due  to  the  variation  of  skin  stiffness 
with  temperature  and  time,   the  Young's  modulus  of  SMA  wire 
is  denoted  by  Ejt).  As  in  Chapter  5,  Equation  (5.4)  could 
be  linearized  according  to  diagram  in  Rogers   [1]   giving  the 
relation  between  the  Young's  modulus  and  temperature.  Since 
the  values  E^{t)  are  different  at  activation  and  deacti- 
vation due  to  hysteresis  effect   (Duerig  et  al .    [2]  p. 98), 
one  can  write : 
At  heating: 

E3(At)    =  3.5  +  0.25    (AT/At)    At  do'  psi)  (7.1) 

E^At)    =  24.13  +  3.10    (AT/At)   At      (Gpa)  T=  °C 

where,   aT  =  32°F  ,  17.78°C 
At  cooling  after  activation: 

E3(At)    =  11.5   -   0.34    (AT/At)    At        (10«  psi)  (7.2) 
E3(At)    =  79.29  -  4.22    (AT/At)   At  (Gpa) 
where,   aT  =  23.5°F,  13.07°C 

Eg  (At)   =  SMA  wire  Young's  modulus 
aT  =  temperature  range 
At  =  time  step  (sec) 
AT/At  a  f(I,V,G)      (Chapter  2) 
1=  electric  current 
V=  electric  Voltage 
G=  heat  flux 
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We  can  approach  the  problem  using  energy  variation  to 
derive  the  equation  of  motion  of  sandwich  panel.  Only  the 
strain  energy  of  the  SMA  hybrid  skins   (upper  and  lower)  and 
the  transverse  shear  energy  of  the  core  are  considered.  We 
ignore  flexural  rigidities  of  the  face  skins  because  their 
thickness  is  very  small  compared  to  the  core  thickness .  The 
recent  research  on  nonlinear  analysis  of  composites  by 
Whitney  and  Leissa   [27] ,  Bai   [28] ,  and  a  collective  work  on 
geometrically  nonlinear  analysis  by  Chia   [29] ,   can  serve  as 
good  models  for  deriving  the  equation  of  motion.  The  effects 
of  skin  stiffness  variation  at  activation/deactivation  will 
also  be  considered. 

The  displacement  fields  are  assumed  in  the  following 

form  (from  a  2-D  plate  approach) : 

Non-activation : 

u(x,y,z,t)  =  Uo(x,y,t)  +  z  i^x  (7.3) 
v(x,y,z,t)  =  Vo(x,y,t)  +z^y  (74) 
w(x,y,t)    =  w(x,  ',t)  (7  5) 

These  assumptions  lead  to  the  conventional  results  found  in 

such  previously  collected  works  as  Whitney  [30] . 

Activation  mode:        Upper  skin  deforms  by  heating  and 

recovery  stresses.   If  the  core  is  an 
insulator  to  heat,   the  extra  terms 
appearing  in   (7.6),  (7.7)   are  small. 

u(x,y,z,t)  =  Uo(x,y,t)  +  z  i^^  -  f(x,t)  (7.6) 
v(x,y,z,t)  =  Vo(x,y,t)  +  z  i^y  -  g(y,t)  (7.7) 
w(x,y,t)    =  w(x,y,t) 
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Where  f  and  g  represent  the  shrinkage  of  upper  skin  at 
activation  (only  at  cantilever  condition) ,  the  compression 
induced  by  SMA  may  be  ignored  if  the  base  structure  is  stiff 
enough,   therby  eliminate  f(x,t)   and  g{y,t)   at  simplv- 
supported  boundary  condition  for  u=0 .   From  Figure  7.1, 
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Figure  7.1.   Cantilever  panel--  extreme  case  of  existing 
f  (x,t)  and  g(y,t)  due  to  upper  skin  activation 


We  introduce  the  core  shear  rotation  due  to  upper  skin 
activation   (cantilever  beam) . 


fix,  t)    -  d 


(7.8) 


where 

Uo(x,y,t)  =  panel  neutral  plane  displacement  in  x-direction. 
Vo(x,y,t)   =  panel  neutral  plane  displacement  in  y-direction. 
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w(x,y,t)     =  panel  neutral  plane  deflection  in  z-direction. 
\l/^{x,y,t)   =  panel  neutral  plane  rotation  about  y-axis. 
l/'y(x,y,t)   =  panel  neutral  plane  rotation  about  x-axis, 
z      =  local  normal  direction  to  neutral  plane  of  panel, 
f (x,t)=  deformation  of  upper  skin  at  activation  in  x-axis. 
g(y,t)=  deformation  of  upper  skin  at  activation  in  y-axis. 
^sx  =  core  shear  rotation  due  to  upper  skin  activation 

in  X  -  direction. 
=  core  shear  rotation  due  to  upper  skin  activation 

in  y  -  direction, 
dx    =        thermal  expansion  in  x-direction. 
dy    =        thermal  expansion  in  y-direction. 

=  skin  thermal  expansion  coefficient  in  x-direction. 
Q!y  =  skin  thermal  expansion  coefficient  in  y-direction. 
To     =        ambient  air  temperature. 

A  diagram  of  the  sandwich  panel,   shown  in  Figure  7.2, 
gives  us  a  clearer  description  of  the  panel's  longitudinal 
displacements  and  rotations.  A  long  beam  with  cantilever  or 
simply- supported  boundary  conditions  under  harmonic  load 
will  produce  large  amplitudes  of  vibration. 

Figure  7.2  only  shows  the  plane  view  of  the  sandwich 
panel  from  the  y-direction;  v„  can  not  be  shown  in  the 
diagram.  Additionally,   dw/dy  can  not  be  shown  except  viewing 
from  the  x-direction,  and  i^x,        are  briefly  denoted  as  il^ 

T  x,y 


Figure  7.2     Large  deformation  of  sandwich  panel 


The  strains  of  skins  can  be  expressed  as 

^  du,  _  hd\P^  ^  1  _  df(x,  t) 
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dx       2  dx       2  ~5x 


(7.10) 


"  2"^  2^W 


dg(y,  t) 
5y 


(7.11) 
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(7.12) 
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(7.13) 


ay        ax       2  dy 


dx  ~3xTy 


(7.14) 


The  nonlinear  parts  of  strains  in   (7.9)   to   (7.14)   are  owing 
to  the  nonlinear  part  of  Lagrangean  strain  tensor  in  Fung 
[31]  : 


^  1  r  du       dv      ,  dudu   ,    dvdv  ^  dwdw,. 


2    ^      Tx        dxdy      dxdy  dxdy' 

3ince  the  second  order  terms  like  (AH)^  .  et( 

ax 

are  small  compared  with  {^)^    ,we  ignore  them. 


The  core  transverse  shear  strains  are 
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Txz  =       -  1/'..  +  (7.15) 


Ty.  =       -        +  4^  (7.16) 


According  to  the  asymmetrical  skin  laminated  sandwich  panel 
at  activation   (upper  skin  stiffness  change)   and  compared  to 
[28] ,   we  find 
skin  strains  expressed 

te"]"  =  [  e",     e"y    y"  ] 

in   (7.9)   to   (7.14)  as 


U']'  -  [  ei     e'y  ] 


Neutral -plane  strain   (which  causes  bending  by  coupling 
stiffness) ,  and  sandwich  panel  curvature  k  are  expressed  as 

[    ,]T  ^    [     d^._di^s.  d^y_Q^Psy  d^yj^.  , 

~Sx      dx        dy      dy       ~Tx  Ty 

Also  according  to  equations  (5.3)  to  (5.11)  in  Chapter  5, 
total  strain  energy  of  the  sandwich  panel  is  expressed  as 


V=  1 


a  b 


12  -"13 


'^1    -^2  -^3 


^31    -^32  -^33 


[  e "]  dx  dy 


a  b 


\l    -^12  Al3 


^31    -^32  -^33 


[e^]  dx  dy 


a  b 


[e    «]  ^ 


dx  dy 


a  b 


h 
■2 


^13  0 

0  G. 


23 


Yxz 
Vyz 


dz  dy  dx 


(7 


Aij  ,         ,  B*  ,  D*  are  the  laminated  plate  stiffness 
properties.  The  superscripts  u  ,   1  ,   c  indicate  upper, 
lower,   core.  The  superscript  *  means  skin  is  activated, 
where 

  N  _ 

^ij  '    E  ^ij  N=upper  skin  layer  number 

♦  :  skin  at  activation 

Aij  :    52  Q".  iz^-z^_^)        N=  lower  skin  layer  number 
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Bii  -613 
B21   .  . 


Dil  Dia 
D21    .  .    .  . 


G,3  : 
G23  : 

Oil  • 


core  transeverse  shear  stiffness  at  xz  plane 
core  transeverse  shear  stiffness  at  yz  plane 
.  .  .etc.  are  function  of      defined  as  eg.   (1)  ,  (2) 


We  could  also  express  the  equation  in  the  following  equation 
(7.18) 


V  =  ± 


•  y 

tr" 


dy  dx 


1 
"2 


]  • 


y 


dy  dx 


+  1 
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a  b 


[  ] 


3  (^^-^«^) 


dy 
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dy  dx 
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a  b 


^  1 


[         gy  ] 


Yyz 


dy  dx 


(7.18) 


The  kinetic  energy  form  could  be  derived  as  follows, 
from   (7.5) ,  (7.6) ,  (7.7) 


u..  = 


u,-^rP^-f(x,  t) 


=  Vo-^^y-9(y,  t) 
w  =  w{x,y,  t) 


a  b 


Poh^  {  ul  +       +       }  dy  dx 


(7.19) 


a  b 


-J  J  J  Po-^i  (  li?  +       +       )  dy  dx 


(7.20) 


a  b 


IJIC    _  1 


pji  {  Uo  +  Vo  +       )  dy  dx 


(7.21) 
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where,  is  SMA  hybrid  laminate  skin  density  (Ib/in^) 

where,  is  sandwich  core  density  (Ib/in^) 

Therefore,  total  kinetic  energy  of  sandwich  panel  is: 
T  =  T"  +        +  T*"  (7.22) 

External  work  can  be  expressed  due  to  a  transverse  uniform 
loading  qo,  as 

a  b 

^  =         [    [  ^0  ^  (7.23) 


Thus,  employing  the  variational  principle  minimizing  the 
total  energy,  we  have  the  following 

U  =  -T  +  V  +  W 

5U  =  -5T  +  6V  +  6W  =  0  (7.24) 

From  Appendix  A  (A. 25) , (A. 26) , (A. 27) , (A. 28) ,  we  assume  that 
under  simply- supported  or  cantilever  boundary  conditions, 
beams  are  under  cylindrical  bending  (1-D  static  condition) . 
Because  we  ignore  the  coupling  stiffness  effect  of  upper, 
and  lower  skins  at  activation,   and  there  are  no  nonlinear 
part  containing  N,"'^  in  the   (A. 28),  we  have  the  following 
governing  equations 
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(7.26) 


The  above  equations  of  motion  are  able  to  lead  to  the  same 
conclusion  as  Whitney  [30]   if  we  neglect   (7.25)    ,    (7.26)  and 
preserve  the  bending  equations   (7 . 27)  ,  (7 . 28)  .  The  \P^^  varies 
in  different  materials  and  is  often  negligible. 

7.2.2     Simply- Supported  Sandwich  Beam 

An  example  case  of  sandwich  beam  under  uniform  loading 
condition   (static)   has  been  solved  by  Whitney  [30] 


=   21^  [        -  2LX^  -  L3  .  ^  (L-x)  ]  (7.30) 


-13^ 
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where  qo  is  the  uniform  load  on  a  sandwich  beam. 
7.2.3     Cantilever  Sandwich  Beam  (static) 
Case  1 : 

A  cantilever  beam,  under  uniform  loading,  q^,   can  be 
solved  by  using  (7.27),  (7.28)   equations  along  with  the 
cantilever  boundary  conditions. 

=  .^2^  .goL^x 

2D,,      -2D^  ^'^.^^) 

detailed  procedures  for  getting  solutions  can  be  found  in 
Appendix  B  . 

Case  2  :  .  f 

A  cantilever  beam,  under  concentrated  load  P  at  the 
tip,   can  be  solved  by  (7.27),  and  the  slightly  different 
Equation   (7.28)   along  with  the  boundary  conditions.  See 
Appendix  B  for  details. 

^    Px^   _  PLx 

"  2D^  "  (7.33) 

^  ^   PLx^  _  Px^   ^  Px 

2D,,      6D„      G,3h  (7.34) 
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7.3    Vibration  Analysis  of  Sandwich  Beam 

The  applications  of  composite  materials  have  been 
spread  from  the  aerospace  engineering  in  the  early  i960  to 
the  automobile  industry,  marine  engineering,   sporting  goods, 
and  some  other  commercial  applications.  We  have  found  many 
research  papers  in  composite  structures,  mostly  related  to 
laminated  beams  and  plates  or  shells  but  very  few  of  them 
are  associated  with  the  studies  of  sandwich  structures. 

Journal  papers  surveyed  by  Kapania  and  Rakesh  [32] , 
[33] ,  have  shown  evidence  of  applying  Classical  Laminate 
Theory  (CLT)   on  thin  plate  analysis  at  early  stages  of 
applications  on  composites.  The  classical  theory,  which 
ignores  the  effect  of  transverse  shear  deformation,  may  not 
exactly  describe  the  kinematics  of  relatively  thick  plates, 
however.   It  can  underpredict  the  deflections  and  overpredict 
the  natural  frequencies  of  thick  plates.  This  is  especially 
true  when  transverse  shear  moduli   (interlaminar  shear 
moduli)  of  laminated  composite  beams  and  plates  are  very  low 
as  compared  to  in-plane  tensile  moduli. 

Therefore,   shear  efferts  are  crucial  in  composite 
laminates.  Following  the  classical  laminate  theory,  two 
categories  in  the  shear  deformation  theory  have  been 
developed.  One  is  the  first -order  theory  and  the  other  is 
the  higher-order  theory.   The  first -order  theory,   the  well- 
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known  Reissner-Mindlin  theory,  was  developed  by  Reissner 
[34]   and  Mindlin   [35] .  More  extensive  work  by  Yang,  Norris, 
and  Stavsky  [36]   and  by  Whitney,  and  Pagano   [37] ,  following 
the  earlier  Mindlin  theory,  has  worked  well  for  laminated 
beams  and  plates  analysis  in  predicting  deflections,  natural 
frequencies,   and  buckling  loads.   The  higher-order  theory,  a 
more  sophisticated  idea,   considers  through-the-thickness 
stress  response,   such  as  cr^,    t^^,   in  regions  of 
discontinuities  at  the  boundaries.  Thus,   the  higher-order 
theory  gives  more  accurate  results  than  CLT  or  Mindlin' s 
theory.  However,  the  theory  requires  much  higher  computing 
costs  and  this  made  it  impractical.  To  overcome  this 
limitation,  Reddy  [38]   has  made  some  progress  to  make  this 
theory  more  applicable. 

The  survey  work  in  [32] ,    [3  3]  also  mentioned  numerous 
other  research  on  vibration  of  composite  beams  and  platt^s 
We  found  references  in  section  7.1  referred  to  this  special 
topic,   in  particular  a  paper  on  the  dynamic  behavior  of 
composites  and  sandwich  platps  by  Bert    [3  9] .   The  papers 
pertained  to  symmetrically  laminated  platPs  and  linear 
vibration  problems  using  the  first-order  or  the  higher-order 
approach.  The  solution  methodologies  involve  the  analytical 
closed  form  solution  as  well  as  the  approximate  numerical 
method  by  Galerkin  and  Rayleigh-Ritz ,   and  numerical 
techniques  such  as  FEM  method. 
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In  recent  years,   there  are  considerable  amounts  of 
research  on  the  unsvmmetric  laminated  composites   (with  non- 
zero B  matrix).  The  reasons  of  this  particular  research  are: 
1)  providing  a  built-in  self-damping  mechanism  when 
subjected  to  dynamic  excitations,   2)   the  results  of 
fabrication  error  in  orienting  different  plies  of  a 
structure  component,   3)   for  aeroelastic  tailoring  to  control 
mode  shapes  and  frequencies  of  lifting  surfaces  or  4)  to 
analyze  the  dynamic  characteristics  of  the  unparalleled 
activations  of  smart  materials  embedded  in  composites  or 
sandwich  structures,   such  as  the  sandwich  beam,  made  of  SMA 
hybrid  composite  skins  and  soft  core.  The  upper  skin  is 
stiffer  at  full  activation  than  the  lower  skin,  indicating 
that         *  0.  Some  analytical  approaches  to  unsymmetric 
laminates  were  presented  by  Whitney   [40]   and  vibration 
analysis  in  unsymmetric  composite  plates  by  Jensen   [41]  have 
laid  a  solid  foundation  in  solving  the  two  governing 
equations  as  shown  in   (A. 27),    (A. 28)   in  Appendix  A. 

Selection  of  the  first-order  shear  deformation  theory 
as  the  assumption  of  sandwich  panel's  displacement  field, 
described  in  section  7.1,   is  appropriate  because  of  the 
panel's  fairly  soft  core  and  the  fact  that  the  thin  SMA 
hybrid  laminated  skins  have  good  bonding  force  with  the 
core.   The  higher-order  theory  may  not  fit  for  unprecedented 
thermal  activation  and  combined  core  deformations. 
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A  long  and  flexible  beam  or  a  sandwich  beam  can  incur 
large-amplitude  vibration,  which  leads  to  nonlinear 
vibration  problems,   so  its  nonlinear  governing  equations 
require  special  technique  to  obtain  fairly  good  results.  The 
researches  in  this  particular  area  can  be  found  in  Kapania  & 
Raciti   [42] .  They  applied  the  FEM  technique  to  obtain  the 
nonlinear  governing  equations,   then  use  the  multiple-scale 
method  to  solve  the  nonlinear  ordinary  differential 
equation.   In  our  case,  however,   the  governing  equations  of 
sandwich  beam  have  already  been  obtained  in  section  7.1. 

Reviewing  the  work  by  Bhimaraddi   [43]    [44]   and  Eslami  & 
Kandil    [45]    [46]   has  helped  the  author  to  solve  a  group  of 
nonlinear  partial  differential  equations.  The  method  assigns 
two  known  admissible  functinn.g  for  w,        which  satisfy  the 
structure  boundary  conditions.  Then  applying  Galerkin  mfithnri 
transforms  the  nonlinear  governing  equations  into  a  non- 
linear ordinary  differential  equation  which  take  the  form  as 
s  +  a^s  +  a2s'  +  ajS^  =  q.  The  final  step  is  to  use  multiple- 
scale  technique  to  solve  this  nonlinear  ordinary  differen- 
tial equation  for  free  vibration  analysis. 

The  latest  research  shows  that  the  study  of  vibrations 
are  mostly  related  to  the  analysis  of  composite  beams, 
plates,   shells  of  solid  laminates.  There  are  very  few 
reports  dealing  with  sandwich  structures,   especially  in  the 
study  of  sandwich  beams  with  embedded  smart  materials.  The 
nonlinear  vibration  problems  of  the  SMA  hybrid  sandwich  beam 
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will  be  studied  for  different  boundary  conditions  in  the 
following  subsections. 

7.3.1     Solution  of  Simplv-Supported  Beam 

The  nonlinearity  induced  by  the  unsymmetric  laminate  of 
the  sandwich  skins  is  expressed  in  appendix  A,    (A. 18), (A. 19) 
and   (A. 22).  The  two  simultaneous  governing  equations   (A. 27), 
(A.  28)    (uo  =  0,  preserving  the  d / dx  {U^dw/ dx)  +  d/dx{N^^dw/dx) 
in   (A.28)),   with  the  substitution  of   (A. 22)   to   (A. 15), 
become 


dw 
T3z 


)  +  A 


+  B. 


dw  d^w 


(7.35) 
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dx" 


(7.36) 


two  nonlinear  partial  differential  equations,  where 


Nx*"     Applied  in-plane  load  due  to  SMA  wire 
recovery  stress  in  upper  skin  (lb/in) 
Nx*^     Applied  in-plane  load  due  to  SMA  wire 


recovery  stress  in  lower  skin  (lb/in) 


I 


Bending  rotational  inertia  (usually  the 
term  is  very  small) 


Equation (7 . 36 )   could  be  expressed  as 
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"THT  ~   T^-nc; — ^  "   n  ^-^      [Non- linear  terms) 

Ox         dx^       G^,h  G^^h  G^^h      dx"  (7  36a) 


Boundary  condition  (simply- supported) 
x=0 

w  =  0  {7.37a) 
a^^/dx  =  0  (7.37b) 

x=L 

w  =  0  (7.37c) 
d^PJdyi  =  0  (7.37d) 
We  could  assume  the  admissible  functions  as 

w  =       S(t)    sin(Q;x)  (7.38) 
lAx  =   (-S(t)/L)   cos(ax)  (7.39) 
a  =  n7r/L 
n  =  1,2,3, .. . 
S(t)   =  flexural  amplitude 
Equations   (7.38)   and   (7.39)   satisfy  boundary  conditions 
described  in  Equations   (7.37a),    (7.37b),    (7.37c),  and 
(7.37d).  Take  the  derivative  of   (7.35)   with  respect  to  x 
then  add   (7.36).   Insert   (7.36a)   in  D,,a/ax  ( a^^/ax)  of 
Equation   (7.35),  we  obtain 
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(7.40) 


D  JL{         +   go  ^  2p,h,+pj2„  _  (n',"^n:')  a^ur, 

Note  that  the  shear  effect  and  the  in-plane  load  effect  are 
included  in  D^^aVax^  term.   Ignore  the  nonlinear  part  in 
DiiaVax^   (...).     Then  substitute   (7.38),    (7.39)    in  (7.40) 
and  apply  Galerkin  method  for   (7.40)  with  weighing  function 
sin(Q!x)  .  The  Equation  for  the  simply- supported  vibration 
beam  (Figure  7.3)  becomes: 


«-  L 


Of  =  nvr/L 

Figure  7.3.     The  simply- supported  beam 


J„   <-s[  ^11  g  ^    )a.'^{K^^K')oL^  jsin^ax 

+  B^^  (a^sin^Qfx-a^cos^axsinofx  ) 
-1  (A,,+A,,)  s^a^cos^axsin^ax  -  |  (A,,-A,,)  ^OL's^Ln^ax  (7.41) 

+  "2  (Aj^-A^J  -^a^cos^axsinofx  >  dx     +       goSinax  dx 
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-s[  D,,(l  +  _f^)a4  +  {K"^N:')a'  ]i:  + 
s=(|c.3B,,  -  )-s3(i^U^^.A,,)   )  (7.41a) 


Eliminating  the  inertia  term,  we  obtain  a  static  problem  of 
external  uniform  load  qo  over  span  of  L.  For  free  vibra- 
tional analysis,  we  delete  q^,   and  (7.41a)   could  be  simply 
expressed  as 


a^s  +  ajS  +  ajS^  +  348^  =  0  (7.42) 

where 


a,  =  {2p,h,+p^h)  (l  +  _JiLa2)   +   la  (7_42a) 


^""-^^'-W^^^' "  (i^;"-i^;^)«^  (7.42b) 


33  =  (I^{A,,-A^^)   -  |B,,cy3)    2  (7.42c) 


^4  =  —iA.^A,,)   4  (7.42d) 


(7.42)   can  also  be  expressed  as 


(7.43) 
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XT*  ^    ht*  ^ 

(2p„h,+p^)  (l  +  ^a^)  + 

^13-" 


Equation   (7.43)   is  Duffing  equation  without  the  damping 
term.  The  current  methods  of  treating  (7.43)   include  Newton- 
Raphson  approach  like  Singh  &  Rao   [47] ,   the  finite  element 
in  time   (FET)   method  of  Bauchan  &  Hong   [48] ,   the  finite 
element  method  of  Hou  &  Yuan   [4  9]   and  the  multiple  .gralP 
technique  in  Bhimaraddi    [43]  ,    [44]  . 

The  cubic  nonlinearity  described  in  (a^/ajs^  is  a 
stiffening  factor  as  the  vibration  amplitude  becomes  large 
enough,  or  this  term  can  also  be  used  to  describe  a 
softening  factor  if   (a,/a,)   is  negative.  This  phenomenon  can 
be  found  in  Cartwell   [50] ,  and  Schmidt  &  Tondl   [51] ,   all  of 
whom  investigate  some  parametrically  excited  systems  or 
coupled  nonlinear  vibration  systems. 

•  3  •  1  •  1 — Case  study  of  linear  vibration 

The  nonlinear  vibration  response  for  a  nonactivated  SMA 
sandwich  beam  or  its  two-faces-activated  case  is  determined 
by   (7.43)   but  with  a3  =  0   ,   since  A,,  =  A,,,   B,,  =  0    (the  case 
is  equivalent  to  a  symmetric  condition).  Therefore,  (7.43) 
becomes 

s  +  w^s  +   (a,/a,)s^  =  0  (7.44) 
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If  we  use  Whittney's  equation   (10.69)    [30,  p.  280]   for  the 
simply- supported  boundary  condition,  we  get  the  linear  eigen 
value  characteristic  equations  defined  in  (7.45),    (7.46)  are 

-iD,,I^*G,,h)  A„  -  G,,hni  B„  =  -ml  A„  .  (7.45) 
-Giai^^  An  -  ^13^^  B„  =  -p^l  B„  (7.46) 

which  is  equivalent  to  the  solution  of   (7.44)  without  the  s' 
term,   in  other  words,   to  the  solution  of  small  amplitude  A. 
For  Mode  1  : 

n=l  ,  wvf  =6.7% 

=  2264.4  Ibf  in  =  2264.4x12  slug  in^  /sec^ 
L      =  12  in 

=  600  Ibf  /in^     =  600x12  slug  /  in  sec^ 
h      =  0.2  in 

I       =  2poh,  x(  (2h,+h) /2)2  =  2.6x10-^    (lb,  in^) 

=  2.6x10-^  /32.2    (slug  in^) 
for  wvf =6. 7%,   1  SMA  hybrid  layer  at  both  top  &  bottom 
face  of  the  sandwich  beam 

p       =  2xpohi  +  p^xh  =  2.675  X  lO'^    (lb„  /in^) 

=  2.675  X  10-2  /32.2  (slug/in^) 
G,3  is  shear  modulus  of   (  XUS  462]fl.nn  )  which  is  blue 


=  0 


(7.47) 


closed  cell  extruded  polystyrene  foam,  produced  by  Dow 
Chemical  Company. 

Dimensional  compatibility  is  very  crucial  for  dimen- 
sional characteristic  equation  in  eigen  value  calculation, 
from   (7.47)   we  get  o}^  =  817.8  rad/sec.  We  could  also  use 
Whittney  [30]   eq.  (10 . 72) ,  (10 . 73)   to  get  the  eigen  value,  for 
1st  mode   (  n  =  1  ) 


CJi  -   


0}^  = 


\2' 


2264.4  X  32.2  X  12     _  ^239.57  rad 


2.675x10-3  sec 


Where,  co,  =  1239.57  x  0.6599    =  817.99    rad/sec  (both  face 
activated  condition  without  in-plane  forces)   comparing  the 
above  approaches  to  the  linear  solution  of   (7.44)   at  null 
ampli-  tude  to  thickness  ratio   (  A/h  =  0) ,   namely  818 .4 
rad/sec,   from  which  we  conclude  that  the  Galerkin  method  is 
accurate  enough.  This  solution  of  nonlinear  modal  analysis 
of   (7.43)  ,  (7.44)   is  based  on  Bhimaraddi    [43]   eq.(18)   cj^,  = 
aj,Il±b,A!±b,A!±b3A^]   +  0(e«)   et  al .  All  the  procedure  is  coded 
in  Gencom2  for  even  further  analysis  of  damped  and  forced 
dynamic  response. 
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7.3.2     Damped  System  under  Harmonic  Excitation 

The  nonlinear  vibration  model,   illustrated  in  the 
previous  section,   is  merely  for  the  study  of  undamped  free 
vibration.   It  has  inherent  features  of  increased  natural 
frequencies  with  respect  to  higher  vibrational  amplitude 
(A/h)  under  the  prescribed  condition.  Our  primary  goal,  in 
this  section,   is  to  excite  the  known  nonlinear  model  with  a 
harmonic  load  to  find  the  dynamic  response  in  either  the 
time  domain  or  the  frequency  domain. 

The  physical  model  is  a  sandwich  beam  of  composite 
faces  with  SMA  wires  embedded,   and  the  sandwich  core  is  a 
good  damping  materials.  All  the  properties  of  the  model  were 
mentioned  in  the  previous  sections.  This  damped  system, 
however,   is  not  a  viscously  damped  but  a  hysteretically 
damped  system.  The  rationale  is  based  on  the  fact  that  the 
composite  material  with  polymer  matrix  follows  the  elastic- 
viscoelastic  correspondence  principle   (its  moduli  can  be 
expressed  as  complex  moduli)  .   In  view  of   [22]   and   [52]  ,  the 
governing  equation  is   (7.42)  with  a  harmonic  forcing  term 
added  on  the  right  hand  side.  The  equation  becomes 

s  +  co^s  +  dzS^  +  djS^  =  eq'   cosf2t  (7.48) 

where        ,   ds,   d^,  were  defined  in  equation   (7.43a)  and 
(7.42a  -  7.42d),    e  is  small  ,   and  q'=  4qo/(7ra,h)    (#/slug  in) 
as  indicated  in   (7.41a).   n  is  forcing  frequency   (rad/sec) ,  t 


is  time  step,  and  s  is  nondimensional  amplitude  (A/h) .  Since 
1  is  very  small  and  L  is  large  for  long  beam,  (7.42a)  can  be 
simplified  as 


a,  =  (2Poh,+p^)  (l  +  -£lLa2)  {7.48a) 


We  then  apply  the  correspondence  principle  from  [53]  and 
[54], 

Gi3  =  Gi3  +  i  Vc  Gi3  (7.48a) 
Dii  =         +  i  r/b  (7.48b) 
(All-All)*  =  (All-All)    +  i  7?e  (All-All)  (7.48c) 
Bii*=Bii  +  i  Vd  Bii  (7.48d) 
(Aii+Aii)*  =(Aii+Aii)    +  i  77a(Aii+Aii)  (7.48e) 
where  Vb  >   Va  '   Vd  >   Ve  are  the  loss  factors  for  the 
composite  laminate   [54] ,   i  =  /-I.   the  values  of  Vb  >        >  Vd 
,   rj^  are  approximately  .  2  ??„     [  t?„  (loss  factor  of  epoxy 
matrix)   =  .0151    ,  core  rj^  is  approximately.    .  12 .  Therefore, 
the  known  values  of  loss  factor  can  be  substituted  in  the 
following  expressions, 


2  ^13-^ 


(7.49) 


d2     =     as  /  ai  (7.50) 

^2  =  -  |BiW)   I  (7.50a) 

dj     =     aj  /  ai  (7.51) 
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^3  = 


3a*L  ,T 


16 


2 
L 


(7.51a) 


From  (7.48a),  the  real  part  and  imaginary  part  can  be 
expressed  as 


From  (7.49),  the  real  part  and  imaginary  part  of  numerator 
can  be  expressed  as 


VGi3+77eGi3)  h 


with 

03^  =  e^  +  i  eg 

es  =   {e^e^  +  e^ej   /  (e^^  +  e^^) 
eg  =   (e^e^  -  e^e^)   /  (e^^  +  e^^) 
from   (7.50a),    (7.51a),   we  get 
^2  =  +  i 

aj  =  a3i  +  i  a32 


i 
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a. 


'22 


a 


■31 


a 


32 


^2    =     (^21^1    +    32262)7(61^    +    62^)  + 

i   (e^a22  -  62321)/ (61'  +  62')   =      e,  +  i  eg 
da  =  (ajiSi  +  33262)7(61^  +  €2')  + 

i   (6ia32  -  63331)  7  (61'  +  62^)   =      69  +  i  eio 

Therefore,   the  nonlinear  ordinary  differential  equation 
(7.48),  with  its  amplitude  s,     follows  the  correspondence 
principle  also  can  be  expressed  as  below  (see   [55]  )  . 

s     =    X  +  i  Y 

We  substitute  this  expression  into  eq  (7.48) 

(X  +  iY)   +   (e^  +  iej  (X  +  iY)    +   (e,  +  ieg)  (X  +  iY)  ^  + 
(eg  +  ieio)  (X  +  iY)3  =  e  q'  cosS2t 


(7.54) 
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Separation  of  the  real  and  imaginary  parts  in  equation 
(7.53)  gives 
real  part 

X  +  ejX  -  e,Y  +  e,{X^-Y^)    -  2esXY  +  e^  (X^  -  3XY^)  + 

e,o(Y^  -  3X^Y)     =(  e  e,q"   )  /  (e.^+e^^)  cosflt  (7.54a) 

imaginary  part 

Y  +  e^X  +  e^Y  +  e8(x2-Y^)    +  2e7XY  -  eg  (Y^  -  3YX2)  + 
e,o(X^  -  3Y^X)     =(-   e  e^q"    )  /  (e.^+e^^)  cosQt  (7.54b) 

where  q"   =  4q^  /  (^rh) 

Solving   (7.54a)   and   (7.54b)   system  of  highly  nonlinear  and 
coupled  equations  with  assumed  initial  conditions,  as 

X(0)    =  X(0)    =  Y(0)    =  Y(0)    =  0  (7.55) 
we  obtained  the  calculated  dynamic  response  under  harmonic 
excitation.   However,   there  is  no  analytical  solution  for  the 
system  of  equations  and  it  is  necessary  to  use  numerical 
method  to  find  solutions. 

Using  the  higher  order  Runge-Kutta  approach  to  obtain 
the  time  response  of  the  physical  model,  we  insert  the 
prescribed  SMA  sandwich  beam  model  in   (7.54a),    (7.54b)  and 
(7.55)   then  apply  ode4  5    (fixed  time  step)   of  MATLAB  to  get 
the  response  plot  as  shown  from  Figures  7.11  to  7.13.  The 
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real  part  X  is  the  physical  nondimensional  amplitude -time 
response  at  a  specific  excited  frequency.  The  frequency 
response  could  also  be  obtained  by  FFT  (fast  fourier 
transform)   over  the  calculated  time  response. 


7.3.3     Analysis  of  Nonlinear  Vibrati 


on 


The  analytical  results  for  sub-section  7.3.1  and  7.3.2 
will  be  discussed  in  this  sub-section.  We  will  discuss  the 
free  vibration  of  two  SMA  hybrid  sandwich  beam  models  with 
different  SMA  wire  volume  fraction  in  each  beam  in  the  next 
sub-sections.  Results  of  a  complete  analysis  of  the  forced 
responses  of  the  two  SMA  sandwich  beam  models  are  discussed 
in  the  sub-section  7.3..  3.3  after  the  free  vibration 
analysis . 

7-3.3.1  Free  Vibration  of  Short  Ream 

The  nonlinear  frequencies   (1st  mode)   of  SMA  hybrid 
sandwich  beam  (s.s)   for  a  short  beam  model  are  described  as 
Model  1   :     L  =12" 

width  =1" 

Three  SMA  volume  fractions  6.7%.   10%,   16%  were  calculated  by 
multiple  scale  method,   and  the  diagrams  are  plotted  with 
three  different  core  depths.   In  Figure  7.4     (wvf)  =0.16 
Figure  7.5   (wvf)   =  Q.i    ,   Figure  7.6     (wvf)    =n . nfi7  there  are 
three  groups  of  curves  represent  the  three  beams  with  core 
depths  0^,    0^,   0^  as  typical  analytical  Model  1  cases. 
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Each  group  has  three  curves  where 

Dash  dot  line(  _._._.)    :  nonactivated  SMA  hybrid  beam 

Dashed  line     (  )    :  upper  face  activated  beam 

S22Ad — line     (  )    :  both  faces  activated 

A     :  free  vibration  amplitude  of  the  midpoint  of  beam 
All  the  plots  were  calculated  under  the  assumption  of  no 
recovery  force  effect    (pure  simply  supported  B.C);  the 
results  checked  with  reference   [56] .  The  same  pattern  of 
nonlinear  frequency  curves  is  shown  for  the  cubic  and 
quadratic  terms  stiffening  effects  while  A/h  is  larger  than 
0^.  The  linear  natural  frequencies   (A/h  smaller  than  0.1) 
of  SMA  hybrid  sandwich  beams  for  three  different  core  depths 
in  Figure  7.4  to  Figure  7.6  show  that  higher  volume 
fractions  of  SMA  gave  lower  natural  frequencies  for  Model  1 
sandwich  beams.  But  the  differences  between  activation  and 
non-  activation  are  more  obvious.  Especially,   the  upper-face 
or  lower- face  activated  condition  at  large  A/h  showed 
nonlinear  frequencies  lower  than  nonactivated  condition. 

However,   the  recovery  force   (NJ   does  appear  in  the 
hinged-hinged  boundary  condition.  The  beam  length  in  the  x- 
direction  does  not  change   (Uo  =  0)  under  this  particular  . 
boundary  condition.   The  activated  SMA  wires  tend  to  shrink 
back  to  the  original  austenite  state  but  are  constrained  by 
hinged  end  support,   and  large  in-plane  forces  are  induced  by 
SMA  wires  which  are  embedded  in  the  composite  laminated 
faces   (fiber  glass/epoxy) . 
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Figure  7.7  gives  the  analytical  results  of  Model  1  beam 
with  SMA  volume  fraction  0.067  as  in  Figure  7.6   ,   but  the 
diagram  clearly  shows  the  increase  of  natural  frequencies 
under  the  effect  of  in-plane  forces  produced  by  SMA  recovery 
forces  in  the  same  beams  as  in  Figure  7.6   .  This  is  the 
hinged-hinged  boundary  condition,  and  the  in-plane  forces 
90#/in  are  estimated  under  the  volume  fraction  of  SMA 
representing  20  wires   (0.01  OD)   at  recovery  stress  58  ksi. 
inserted  in  2  plies  of  fiber-glass  prepreg.  Each  wire 
creates  approximately  4 . 5  bounds  of  force . 

7.3.3.2     Free  Vibration  of  Long  Beam 

The  nonlinear  frequencies   (1st  mode)   of  SMA  hybrid 
sandwich  beam  (s.s)   for  a  long  beam  model  described  as 
Model  2    :     L  =120" 

width  =1" 

Two  SMA  volume  fractions  6.7%.   16%  were  calculated  by  the 
multiple  scale  method,   and  the  diagrams  are  plotted  with 
three  different  core  depths.   In  Figure  7.8     (wvf)   =  n  ifi, 
Figure  7.9   (wvf)   =  0.067  ,   Figure  7.10     (wvf)    =  0.067   .  -in- 
plane  force  effect.  There  are  three  groups  of  curves,  in 
each  figure,   representing  the  three  beams  with  core  depths 
0-^/  0^  as  typical  analytical  Model  2  cases.  Each 

group  has  three  curves 

Dash  dot  line(  _._._.)    :  nonactivated  SMA  hybrid  beam 
Dashed  line     (  )    :  upper  face  activated  beam 
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Solid  line     (  )    :  both  faces  activated 

A     :  free  vibration  amplitude  of  the  midpoint  of  beam 
The  inplane  force   (45#/in)   in  the  face  sheet  of  the  sandwich 
beams  of  Figure  7.10  was  assumed  value,   it  is  due  to  the 
possible  slippage  of  wires  if  the  ends  at  boundary  fixing 
devices  are  not  tight  enough. 
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Figure  7.4    Free  vibration  nonlinear  frequency  of  Model  1  ] 
(wvf)    =16%  > 

"I 
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Figure 


7.5,  Free  vibration  of  nonlinear  frequency 
(wvf)   =  10% 


Model 


1 


Figure  7.6.  Free  vibration  of  nonlinear  frequency  Model  1 
(wvf)   =  6.7% 


Figure  7.7.  Free  vibration  of  nonlinear  frequency  Model  1 
(wvf)   =  6.7%     (inplane  forces  effect) 


Figure  7.8    Nonlinear  frequency  of  Model  2  (wvf)=l6% 
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Figure  7.9     Nonlinear  frequency  of  Model  2  (wvf)=6.7% 
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Figure  7.10  Nonlinear  frequency  of  Model  2  (wvf)=6.7% 
(inplane  force  effect) 


134 

In  comparison  to  Figure  7 . 6  and  Figure  7.7  ,   Figure  7 . 9  and 
Figure  7.10  manifest  greater  differences  in  recovery  stress 
(in-plane  force)  effects.  The  design  of  SMA  hybrid  sandwich 
beam  could  utilize  the  said  effects  in  vibration  suppression 
for  short  and  long  beams,  and  adding  the  good  damping  effect 
of  sandwich  core  materials  also  reduces  acoustic  noise  in 
many  applications. 

Due  to  the  fact  that  design  of  structure  components 
seeks  to  obtain  lighter  weight  and  lower  stress  levels,  a 
smaller  SMA  volume  fraction  (but  not  smaller  than  3%)   in  the 
laminate  is  recommended.  The  internal  stresses  induced  in 
laminates,  while  the  wires  are  activated,   can  be  decreased 
if  a  lower  volume  fraction  of  SMA  wires  are  embedded. 

'7 .2  .3  .3 — Dynamic  Response  of  SMA  Sandwich  Beam 

The  dynamic  responses  of  the  simply-supported  beam, 
described  in  the  sub-section  7.3.2,   are  calculated  and 
plotted  in  the  following  Figure  7.11,   7.12,  7.13. 

Each  diagram  has  an  upper  section  and  a  lower  section; 
the  upper  represents  the  time  response,   and  the  lower  is  its 
response  in  the  frequency  domain  under  a  particular  exci- 
tation frequency. 

The  text  "omega"  at  the  left  side  of  each  diagram  is 
the  harmonic  forcing  frequency.  The  y- label  "  A/h  " 
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represents  the  " (deflection  of  middle  point  of  the  beam)  / 
(thickness  of  the  core) . " 

In  Figure  7.11,   the  time  response  of  a  long  beam,  under 
the  harmonic  excitation   (omega=0 . 6366  HZ  =  4  rad/sec)   at  the 
beam  center;  showing  vibrational  amplitudes  of  this  beam's 
center  at  the  nonactivated  (N^*"=N/^=0)   and  the  activated 
(N^*"=N/^=45#/in)   condition.  We  use  4th/5th  order  Runge-Kutta 
-Fehlberg  method  to  calculate  the  time  response  of  Equation 
(7.54a)   and  (7.54b)   from  to=0.0  second  to  tf=6.0  seconds. 
Totally  1024  time  steps  were  taken  and  each  time  interval 
were  fixed.  Use  of  the  Manning  window  on  these  two  steady- 
state  responses,  we  can  apply  fast  fourier  transform  (EFT) 
to  transform  dynamic  responses  from  time  domain  (upper 
diagram)   to  frequency  domain   (lower  diagram) .   The  lower 
diagram  shows  response  amplitudes  of  these  nonactivated  and 
activated  beams,   it  also  show  the  excited  frequency  and  the 
natural  frequency  of  each  condition.  From  the  lower  diagram, 
the  nonactivated  beam's  response  has  two  peaks;  one  is  at 
omega=.6366  HZ  and  the  other  is  at  beam's  natural  frequency 
which  can  be  cross-checked  in  Figure  7.10.  The  activated 
beam's  response  also  possess  two  peaks  where  one  peak  is  at 
excited  frequency  and  the  other  peak  is  at  the  natural 
frequency,  which  is  found  compatible  to  that  of  Figure  7.10. 

In  Figure  7.12,   responses  of  the  same  beam  (model  2), 
under  the  harmonic  excitation   (omega=3.183  HZ  =  20  rad/sec); 
at  the  same  nonactivated/activated  conditions  and  solution 


methods.  Because  the  excited  frequency  is  close  to  the 
nonactivated  beam's  natural  frequency,   the  beat  phenomenon 
happens  for  this  nonactivated  beam.  The  lower  diagram  shows 
that  the  natural  frequency  increases  for  this  nonactivated 
beam  as  the  beam's  vibrational  amplitude  increases.  The 
diagram  also  has  two  major  peaks  and  other  small  peaks  in 
frequency  domain.  This  phenomenon  may  be  resulted  from  its 
nonlinear  behavior  at  near-resonance  occurrence.  Also,  the 
activated  beam's  vibrational  amplitude  has  two  peaks 
represent  the  excited  frequency  and  its  natural  frequency. 

In  Figure  7.13,  the  same  nonactivated  and  activated 
conditions  for  this  long  beam  were  under  near-resonance 
excitation  which  is  close  to  the  natural  frequency  of  this 
activated  beam.  The  solid  line  is  the  time  response  for 
activated  beam,   the  response  also  shows  beat  phenomenon.  The 
two  peaks  of  solid  line,   in  frequency  domain,   represent  both 
the  excited  frequency  and  natural  frequency.  The  natural 
frequency  is  about  6  HZ,  greater  than  5  HZ  depicted  in 
Figure  7.10,  which  is  due  to  its  increasing  vibrational 
amplitude . 

Apparently,  under  this  harmonic  excitation  (omega= 
5.013HZ),   we  can  deactivate  the  activated  beam  by  cooling 
it.  The  beam  transforms  back  to  nonactivated  condition  as 
the  dashed  line  represents  in  Figure  7.13.   This  is  how  we 
can  tune  the  structure  to  shift  its  natural  frequency. 
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Figure  7.11.     Dynamic  response  of  Model  2  (q"=.l) 
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Figure  7.12.     Dynamic  response  of  Model  2  (q"=.i) 
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Figure  7.13.     Dynamic  response  of  Model  2  (q"=.l) 


CHAPTER  8 
CONCLUSION  AND  FUTURE  WORK 

The  current  problems  in  dealing  with  the  vibration 
control  of  structures  are  primarily  attributable  to  two 
concerns.  One  is  to  reduce  the  stress  level  in  structural  or 
mechanical  parts  when  the  vibration  occurs  by  methods  such 
as  damping  to  reduce  the  vibrational  amplitude  to  achieve 
longer  operational  life;  the  second  is  to  control  the  modal 
response  of  structures.  Using  SMA  in  the  composite  structure 
is  part  of  the  solution  to  both  problems. 

8 . 1     Conclusion  of  Current  Study 

There  are  five  major  research  areas  in  this 
dissertation.  They  are: 

•  Heat  conduction,   transient  thermal  responsp. 

We  can  make  a  structure  stiffer  or  more  flexible  by 
adding  or  dissipating  a  small  amount  of  thermal  energy. 
Chapter  1,   Chapter  2,  and  Chapter  4  describe  SMA  properties 
and  the  transient  thermal  response  with  test  comparisons,  and 
analyses  so  that  the  activation/deactivation  process  in 
hybrid  structures  can  be  analyzed. 

•  Stress-strain  relations  of  SMA  hybrid  composites 

We  have  to  deal  with  this  problem  with  care  because  of 
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its  complexity.   In  a  composite  plate,   the  stress-strain 
relations  must  combine  with  the  internal  stresses  induced  by 
the  SMA  wires.  This  is  complicated  by  the  two  distinct 
behaviors  of  SMA  at  the  two  different  material  phases,  hot 
and  cold,  which  has  generated  the  special  term,  recovery 
stress .  due  to  initial  straining  before  embedding  into  the 
structure.  The  author  used  a  simple  modeling  approach  to 
achieve  near-uniform  internal  stress  distribution  in 
composite  laminate   (0°  fiber  alignment  using  SMA  wires  as  a 
first  case) .  We  can  simply  use  classical  laminate  theory 
plus  internal  stress  distribution  analysis  to  estimate  the 
laminate  strength  during  the  activation/deactivation .  This 
approach  is  still  in  the  scope  of  macromechanics . 

•  Making  of  SMA  hybrid  composites 

The  fabrication  of  even  a  simple  SMA  laminate  beam  of 
four  to  eight  layers  is  tedious  work,   especially  in  the 
tooling  for  SMA  wire  pre -straining  and  fixing  which  has  many 
more  procedures  than  that  of  plain  laminates.  The 
description  of  this  is  in  Chapter  3 . 

•  Testing  of  specimens 

This  is  another  time-consuming  process  for  SMA  hybrid 
structures.  The  boundary  conditions  for  the  hybrid 
structures  must  be  compatible  with  theoretical  assumptions. 
The  experimenter  must  avoid  data  errors  and  imperfect 
specimens  which  can  make  the  comparison  impossible.  The 


experiments  described  in  Chapter  4  and  Chapter  6  show 

different  levels  of  complexities. 

•        Dynamic  behavior  of  SMA  hybrid  structures 

The  equation  of  motion  of  a  sandwich  beam  in  general 
has  extended  to  the  first  order  nonlinear  approach.  The 
assumptions  are  based  on  large  deformation  theory  and  the 
hinged  support   (hinged- simply  supported)  boundary  condition 
which  could  be  modeled  conservatively  as  a  beam  or  a  panel 
fastened  to  a  base  structure. 

This  particular  hinged  support  beam  case  can  also  be 
applied  to  a  2-D  sandwich  panel.  We  just  add  more  related 
terms  in  the  2-D  governing  equations,   and  the  admissible 
functions  can  be  obtained  in  like  fashion. 

Basically,  we  used  the  natural  frequency  shifting 
method,  an  open  loop  control  system,   to  suppress  vibrations. 
The  stiffening  effect  is  greater  for  a  sandwich  panel  with  a 
thicker  core   (diagrams  in  Chapter  7)   even  without  the 
recovery  stress  effect.  This  is  even  more  valid  for  softer 
face  sheets. 

The  application  of  SMA  hybrid  structure  with  a 
cantilever  boundary  condition  is  also  possible.  We  can  use 
active  property  tuning   (APT)   to  increase  its  flexural 
stiffness  to  control  the  static  deflection  or  buckling 
strength.  Another  approach  for  modal  control  of  a  rotating 
cantilever  structure  such  as  a  large  rotating  wind  turbine 
could  be  active  strain  energy  tuning   (ASET) . 
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8^_2  Future  Study  of  SMA  Hybrid  Structures 

Suggestions  for  future  study  include: 
Use  in  short  fiber  composites  and  other  materials. 
The  base  composite  material  may  not  have  to  be  uni- 
directional fiber-reinforced  composites.  They  can  be 
short  fiber  (SMC)  panels,  plain  epoxy,   and  plastics. 
Although,   coating  of  SMA  wires  might  be  required  for 
better  bonding  with  base  composite  materials. 
Prediction  of  microstresses  in  SMA  hybrid  compos-i  tea . 
SMA  hybrid  structures  may  have  higher  internal  stresses 
with  higher  wire  volume  fraction,  especially  at  a 
cantilever  boundary  condition  as  depicted  in  last  few 
diagrams  in  the  Chapter  5.  The  o^,   t„  (composite's 
compressive  stress  and  shear  stress)   combined  with  the 
hygrothermal  effect  of  the  laminate  can  generate  very 
complicated  problems  in  micro-crack,   resulting  in  local 
failure  at  the  microstress  level.  The  shear  stress  and 
peeling  stress  at  the  bonding  interfaces  between  the 
face  laminates  and  sandwich  core  may  cause  debonding  at 
the  activation.  Large  amplitude  vibrations  of  a 
sandwich  beam  may  obtain  stronger  bonding  force  by 
using  high  performance  film  adhesives.   In  the  study  of 
fracture  problems,   the  activation  of  SMA  wires  in 
composites  may  result  in  the  closure  of  small  cracks. 
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The  use  of  SMA  wires  in  woven  composite. 
Stitching  method  is  potentially  very  helpful  to 
forestall  delamination  and  improve  dynamic  properties 
simultaneously.   In  these  cases,   the  modeling  of  the 
hybrid  structure's  stiffness,   strength,  microstress  is 
very  critical. 

Other  potential  applications. 

In  modal  analysis,  using  SMA  to  control  low  frequency 
vibrations  of  long  flexible  beams  or  panels  made  of 
hybrid  composites  could  be  very  important . 
Using  SMA  as  both  sensor  and  actuator  in  more  advanced 
closed  loop  vibration  control  devices  for  long  flexible 
structures . 

Using  SMA  in  structures  for  acoustic  noise  control. 
And,   in  configuration  problems,  using  SMA  in  structures 
to  control  the  geometric  shapes. 


APPENDIX  A 
THE  EQUATION  OF  MOTION  OF  SANDWICH  PANEL 


The  variational  principle  of  (7.24)   in  Chapter  7  will 
be  derived  in  details  as  the  following: 


a  b 

dT  = 


j  J  Po^i  5{  ul  +  vi  +       )  dy  dx 


a  b 


+  ^  J  J  Po-^i  5{  uj  +  vj  +       )  dy  dx 


a  b 


+  ^  J  I  Pc^  5  (  lio'  +  v^  +       )  dy  dx  (A.l) 
For  example,   taking  the  upper  skin  velocity  in  x-dir.  we  get 

^     ^      *  a  b 

2  I  I  ^"^^  dy  dx  =         J  J 


a  b 


-•'''J   [    ^^^'^-^'^rdx  ^  ^^^^^ 


Po-h,  J  J  (-g^         \  -  ^  6uJ  dydx 
6ujO}  =6ujt,)  =0 


assumed  that  the 
same  fasion  for 

i 
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Therefore, 

a  b 

=  j  J  ^-P^^i       S".  -PoK  ^1  <5"i  -PoK  -Po^  ^i^^i     (A.  3) 


since, 


2  '  X  ---OS 

....  h 

2  y  ""'0  "2"^^' 


After  several  manupulations,    (A. 3)  is 


a  b 


(  -(2p„h,+p^h)      +  p^hj  )  6v, 
(  -  (2p(,h,  +  p^h)  w'  )  dw 

(  -p,h,u,--^ip^-p^hf  )  af 
(  -p,h,v^-h^ip  -p^h,g  )  6g  ]  dydx 


(A. 4) 


From  Chapter  7,   equation   (7.23),   we  get 


a 

"  "I  J   j  ^^"^  (A.  5) 


From  Chapter  7,  equation  (7.18),  we  also  can  apply  the 
variational  principle  to  get 


b 
a 

a 
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a 
a 


a 


b  a 


ax 


ag; 

-3 —  )  5w  ]  dydx 
dy 


(A. 6) 
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Since  the  activation  term  in  (A. 4)   is  not  accelerating 
quickly  enough,  we  ignore 

f  term  ,  {  f  =  h  il/^J 
g  term  ,  {  g  =  hip^^) 

Therfore,   from  (A. 4),    (A. 5),  and  (A. 6),  we  have  the  line 
integrals  representing  the  boundary  conditions  of  the  panel. 
Then  we  collect  the  terms  having  the  same  variation  terms  as 
6uq  ,   6vo,  and  so  forth,   in  the  double  integrals.  We 
conclude 
6un  term  =  0 


6v„  term  =  0 

^-J3r  ^  ^  '  =  (2Po^.p,h)  v„  (A.8) 


d\l/^  term  =  0 
2       dx        dx  ^yF 
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6\p    term  =  0 


(A. 10) 


term  =  0     ,    {5\p^^=   (i/h)  6f) 
""^-^^  -  (A.ll) 


May  term  =  0     ,    (5ip^^=   (l/h)  6g) 


6w  term  =  0 


^        F        ^  "^^"^^^^^^"^^^^ 

TJc   -  ^^""^^^    "53>^^^a3^^  ^^^^-^^       (A. 13) 


ag;  ag; 


Assuming  1-D  beam  model 
In  the  meantime, 
we  assume  • 

So,  we  only  have  4  equations 
they  become 


,  we  ignore  N/-^  ,  N^"'^  terms. 
(A. 7) , (A. 9) , (A.ll) , (A. 13)  left. 
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dN"  ^  awj 

dx  dx 


=  0  (A. 14) 


(A. 17) 


In  Chapter  7,  equations  (7.9)  -  (7.17),  we  are  able  to 
derive  the  following  expression 

dx'       2  ax^       a3?  ~^JJ^ 

axay     ax^     2  axay     2  T?" 

^  dw  a^K^  _^  av  a^v  , 

Ty- ax^  dxdy 


The  Aij  terms  are  the  upper  skin  stiffness  in  activation, 
and  the  Ai^  terms  are  the  lower  skin  stiffness. 
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dx^       2  dx^       'Jx  dx^ 

axay     ax^     2  "5x5y    2  ax^ 

dw  d^w  ^  avv^  a^w  J 
ax^  dxdy 

^  =  A     (  -  -  +   dw  d'w  , 

dy        "     dxdy      2  dxdy      dxdy  dyT3^ 

+  A,  -  ^i!^2f  +  dwd^  _  a^g 

dy^       2  dy^       dydy-       a^  ^ 


dy'       dxdy      "2  axay      "2  ~a>^ 

^   a^^  a^w  ^   dw  d^w  J 
"5x  ay^     Ty  axay 

dy^        2  dy^       '^dy^'  > 

.A33         .  ^^^0  ,  h  a^^,  ^  ^a^^^  (^-2^) 

ay2       '^xUy      2  'J3c^      2  ~dy^ 
^  dw  d'w  ^   an^  d^w  J 
ay^  dxdy 

In  the  same  fashion,  by  using  (7.11  -  7.17),  we  get 

ax     ax    ""^^^^  ^  ^    -^^^^^  -rf^ 

.  aiAy  a,A,  au    1  a^..^, 

^^^^  2^^^  ^  ^ 
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dy     dy  ^  ""^'^^  "aF  ^    ^"^-57  "aF 


gx  =  G,3hyL  =  G,3h(^,  -        +  ||)  (A. 24) 


Assuming  the    A^^=A^^,  B^j  =  0,  at  a  moderate  activation, 
no  resultant  shear  N^"'^  no  twisting  term  as         in  1-D 
beam,   substituting   (A. 18) , (A. 19) , (A. 22) , (A. 23) , (A. 24)  into 
the  previous  equations  of  motion   (A. 14),    (A. 15),    (A. 16), 

(A. 17).  We  conclude  that  in   (A. 18),   all  A^^  (  )   terms  in 

the  bracket  are  eliminated  for  the  1-D  beam  case.  Same  as  in 

(A.  19),   all  A21  (  )   are  eliminated.  The  SMA  lay-up  in 

laminate  is  considered  the  same  as  the  orthotropic  laminates 
mentioned  in  Chapter  5.  Then,   the  A31  in   (A.  18),   and  the  A31 
in   (A. 19)   are  zeros,   too.  Thus,   Equation   (A. 14)  becomes 


2  A,,{ 


2  dx^ 


) 


(A. 25) 
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Since,  ignoring  the  sandwich  beam  twisting  moment  M^, 
(A. 23)   could  be  eliminated, 

and  ^.  a^^^  ^  Dl^df 

dx  h 


D2i(...)   and  D31  (....)   in  (A.  22)   are  eliminated  for  the  1-D 
case.  Therefore,    (A. 15) , (A. 16) , (A. 17)  becomes 

G,ArP    -rp     ^  =  -  (A  26) 

i3^V.     Vs.  -j^-j^      -j^j^,  IA.26) 


A,  A   1 


dw  d^w 


2  dx^ 


dx^ 


)      G,A\p    -  \P     +  4^) 

13  '     X         r  3x  4  ' 

h  dx"  a3? 


(A. 27) 


The  f (x,t)   and  g(y,t)   are  very  small  and  can  be  treated 
as  perturbed  terms;  additionally,   they  are  different  in 
various  kinds  of  SMA  composites. 


APPENDIX  B 
CANTILEVER  SANDWICH  BEAM 


Case  1 


■■■■■■■lllilllllllllilll 


/ 
/ 

A  cantilever  beam  under  uniform  loading  qo  with 
governing  equations 


(B.l) 


and 


.  ^     ^     aP^     ^°  "  ° 

For  simplicity,  we  define  -  \p^^  by  so  equations 
(B.l) ,  (B.2)  become 


(B.2) 


h  d: 


(B.l) 


and 


=  0 


(B.2) 
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Boundary  conditions  are 
at  x=0 


w=0         ,         rl^^=o  (2  3) 

D,,d\PJdx  =  -  qoLV2  (B.4) 
Giahi^^-x  +  aw/ax)    =  qo  L  (b.5) 


and  at  x=L 


d\PJdx  =  0  (  =  0)  (B.6) 
lAx  +aw/ax  =  0     (  shear  =0  )  (B.7) 

Integration  of  eq.    (B.2)   and  using  boundary  condition  (B.5) 

we  obtain 

Gi3  +  11)=  -go^  ^<Io^  (B.8) 


Substitution  of  eq.    (B.8)   into   (B.l)   leads  to 


D. 


11 


The  second  eq.  of   (B.9)   is  obtained  by  integrating  the  first 
eq.   of   (B.9)   and  from  boundary  condition   (B.6),  we  get 
C2  =  -  qo  I//2.  Then,   substituting  the  second  equation  of 
(B.9)   and  Cj  into   (B.2),   we  obtain 


d^w  _ 

g,Lx  ^ 

2-Dii 

^11 

2D,, 

dw 

-  \^  - 

go^' 

6i^ii 

q,Lx^  ^ 
20,, 

2D,, 

(B.IO) 
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Based  on  the  boundary  condition  (B.3),   the  second  equation 
of   (B.9)  becomes 


2D, 


2  A 


X  +  5 


(B.ll) 


The  term  6  is         (the  core  shear  rotation  due  to 
activation) . 

For  the  deflection  w(x) ,  we  substitute  (B.IO)  and 
(B..11)   into  boundary  condition  (B.7).  This  leads  to 
C3=qoL/Gi3h 
and 

w  =  ^  -S^  M^-^)x^  ^^x  .6x  ( 


Case  2 :  /l^  L  -» 

/I  t 

P 

For  a  cantilever  beam  under  concentrated  load  P  at  the 
tip,   the  governing  equations  are 
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The  boundary  conditions  are 
x=0 

w  =0  ,        =  0 

Dii  d\pjdx  =  -PL 


dxl/Jdx  =  0 

lAx  +  dvi/dx  =  P/Gijh 


(B.15) 
(B.16) 

(B.  17) 
(B.18) 


integrate  (B.14) 

Gi3h   (        +  aw/ax  )   =  Ci  (B.19) 

and  substitute  into   (B.13),   it  becomes 

=  C,        /2D,i  +  C2  X  (B.20) 

which  satisfies  the  boundary  condition  (B.15).  Using  the 
boundary  condition   (B.17),   we  obtain        =  -  C^L/  D^^. 
Therefore,        can  be  expressed  as 

xp^  =(  xV2Dii   -  Lx/Dii   )Ci  (B.21) 

Substituting  (B.21)   into   (B.14)  we  obtain 
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a  r  2  B.22 

=        -  _£L)c  +  c 
ax     'd,,  idJ"-^ 


From  eqs   (B.21),    (B.22)   and  the  boundary  condition  {B.18), 
we  get 

C4  =     P  /  Gijh  {B.23) 

take  the  derivative  of  eq  (B.21)  with  respect  to  x  and  apply 
the  boundary  condition  (B.16)   to  derive 

Ci  =  P  (B.24) 
Finally,  we  substitute  the  boundary  condition  (B.15) 
and  eqs   (B.23),    (B.24)   back  into  eqs   (B.21),    {B.22)   with  the 
integration  of  the  second  eq  of   (B.22),  and  conclude  the 
following: 


_    Px2       PLx  . 

-   2^  "  -D^T  (B.25) 

„  _   PLx^      Px^    ^  Px  . 

-  3^  ~-6D^  "  G:;^  ^  (B.26) 


The  term  d  is  \J/^^  (the  core  shear  rotation  due  to  the 
sandwich  beam  activation) . 
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